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ABSTRACT 

The  objective  of  this  study  was  to  find  optimum  design  features  for  deep 
underground  protective  structures  in  rock.  A  computer  program  was  developed  to 
determine  survival  distance,  optimum  depth,  total  construction  costs,  and  costs  per 
usable  volume  for  protective  structures  in  rock.  The  program  was  applied  to  structures 
with  inner  radii  of  3,  6,  and  9  meters:  concrete  liners  0.5,  1,  2,  and  3  meters  thick 
and  a  steel  liner  2.54  centimeters  thick  were  investigated.  In  the -example,  the  rock 
field  was  assumed  to  be  either  of  sandstone  or  granite.  The  influence  of  the  effective 
longitudinal  seismic  velocity  between  the  detonation  point  and  the  structure  was 
studied.  Survival  distances,  survivabilities,  optimum  depths,  and  pertinent  costs  were 
found  for  true  surface  bursts  of  100  and  1,000  kilotons. 

It  was  found  that  a  concrete  structure  in  sandstone  can  provide  full  protection 
against  a  weapon  yield  no  larger  than  100  kilotons,  and  then  only  if  the  structure  has 
a  sufficient  soil  cover.  Full  protection  against  a  weapon  yield  of  100  kilotons  is  pro¬ 
vided  by  a  structure  with  a  steel  liner  in  granite  irrespective  of  the  composition  and 
stratification  of  the  rock  above.  Full  protection  against  a  weapon  yield  of  1,000  kilotons 
can  be  provided  by  structures  with  steel  liners  in  granite,  at  depths  less  than  300  meters, 
only  when  the  effective  longitudinal  seismic  velocity  between  the  detonation  point  and 
structure  is  2,000  meters  per  second,  or  less.  In  all  the  cases  where  full  protection 
cannot  be  achieved,  the  optimum  depth  it  between  100  and  300  meters. 
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INTRODUCTION 


The  purpose  of  this  report  is  to: 

1.  Obtain,  from  fundamental  solutions  of  mechanics  and  available 
test  results,  data  on  the  protective  ability  of  structures  in 
rock, 

2.  Determine  the  construction  costs,  and 

3.  Establish  the  costs  for  various  levels  of  protection. 

The  protective  ability  of  a  structure  can  be  measured  by  the  distance 
between  the  structure  and  the  point  of  a  nuclear  detonation.  If  the  ac¬ 
curacy  of  the  attacking  weapon  system  is  taken  into  account,  the  pro¬ 
tective  ability  can  be  expressed  by  the  probability  of  survival  (the  sur¬ 
vivability).  Both  the  survival  distance  and  the  survivability  are  depend¬ 
ent  upon: 

1.  Characteristics  of  the  attacking  weapons, 

2.  Geological  environment,  and 

3.  Structural  design. 

These  factors  are  embodied  in  a  computer  program  and  an 
example  configuration  is  used  to  illustrate  the  influence  of  the  princi¬ 
pal  parameters  on  the  survival  distance,  survivability,  and  cost  effec¬ 
tiveness. 

Different  optimum  configurations  exist  for  underground  protec¬ 
tive  structures  because  of  various  weapon  characteristics,  geological 
features,  and  structural  geometries.  The  trend  toward  optimum  con¬ 
ditions  has  to  be  found  to  establish  guidelines  for  site  selection  and 
design  and  to  determine  the  cost  effectiveness  for  underground  protec¬ 
tive  structures  in  rock,  based  upon  military  requirements. 

For  protection,  the  guidelines  for  site  selection  should  contain 
numerical  values  for  the  significance  of  rock  types  and  tectonic  con¬ 
struction  patterns.  Guidelines  for  designing  underground  protective 
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structures  should  contain  the  optimum  depth,  size  of  cross  section,  type 
of  liner  material,  and  liner  thickness  for  the  various  cases.  The  eco¬ 
nomic  evaluation  should  indicate,  for  a  particular  protective  purpose, 
whether  it  is  better  to  have  one  structure  of  a  higher  single-shot  sur¬ 
vivability  and  higher  unit  cost,  or  to  have  several  structures  of  lower 
single- shot  survivability  and  lower  costs.  An  economic  evaluation  also 
might  determine,  to  some  extent,  the  most  economical  direction  for 
further  studies.  No  single  source  containing  solutions  or  information 
on  the  entire  problem  was  found  in  a  survey  of  available  unclassified 
literature  on  the  design  of  protective  structures  in  rock.  However, 
many  books,  reports,  and  papers  together  contained  all  the  parts  of 
the  problem.  The  main  informational  sources  for  this  research  were: 

1.  Reference  1— rock  mechanics, 

2.  Reference  2— statically  caused  stress  distribution  around 
openings,  and 

3.  Reference  3— dynamically  caused  stress  distribution  around 
openings. 

Since  References  2  and  3  were  basic  to  the  information  in  this  re¬ 
port,  both  static  and  dynamic  load  on  underground  protective  structures 
were  considered  so  the  influence  of  the  structural  depth  could  be  studied. 
However,  the  application  of  References  2  and  3  placed  an  important 
limitation  on  this  investigation;  that  is,  that  the  liner  material  and  the 
surrounding  rock  up  to  a  thickness  of  3  cavity  diameters  is  assumed 
to  behave  elastically. 

In  spite  of  this  limitation,  the  computational  procedure  based  on 
the  elastic  behavior  of  material  is  of  theoretical  and  practical  signifi¬ 
cance.  It  is  the  first  step  toward  improved  procedures  that  may  include 
the  plastic  range,  and  it  can  be  used  directly  for  designing  protective 
structures  capable  of  resisting  several  attacks.  The  second  limitation 
placed  on  this  study  is  that  only  the  directly  induced  ground  shock  is 
taken  into  account.  Consequently,  this  report  deals  only  with  pro¬ 
tective  structures  located  at  depths  in  which  only  the  directly  induced 
ground  shock  is  critical. 

Optimum  conditions  are  ascertained  by  seven  major  steps: 

1.  Determination  of  the  stress  distribution  around  the  opening 
in  the  liner  and  in  the  surrounding  rock, 
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2.  Determination  of  the  stress-dependent  strength  of  the  liner 
and  the  surrounding  rock, 

3.  Determination  of  the  survival  distance, 

4.  Determination  of  the  survivability,  based  on  survival  distance 
between  detonation  point  and  structure, 

5.  Determination  of  the  gross  construction  costs  and  costs  per 
usable  volume  unit,  based  on  the  allowable  design  dimensions 
of  the  structures  and  the  construction  procedures, 

6.  Assembly  of  a  computer  program,  and 

7.  Determination  of  optimum  design  conditions  . 

THEORY 

Dstsrminatton  of  Stross  Distribution 

The  load  on  an  underground  protective  structure  consists  of  re¬ 
sidual  forces  in  the  rock  body,  the  weight  of  the  overburden,  and  the 
ground  shock.  Since  they  are  unpredictable,  residual  forces  have  been 
neglected  for  this  feature  study.  However,  residual  forces  may  be  in¬ 
troduced,  as  in  overburden,  into  the  stress  computation  if  they  can  be 
measured  in  a  particular  case. 

The  weight  of  the  overburden  is  to  be  considered  not  only  because 
of  its  contribution  to  the  stress  distribution  around  the  opening  but  also 
because  of  its  effects  on  the  strength  of  rock  and  liner  materials. 


Determination  of  Stresses  Caused  by  Overburden.  To  compute 
the  stress  distribution  caused  by  the  static  load  of  the  overburden 
weight,  Savin's  solution  was  used.  His  solution  is  the  only  one  known 
that  considers  a  plane  distributed  load  on  a  lined  circular  opening  in 
an  elastic  field.  Other  solutions,  as  in  Reference  1,  are  related  to 
a  circular  distributed  load. 

The  equations  of  the  stresses  are  given  in  Appendix  A.  The  in¬ 
fluencing  factors  are  depth  of  the  structure,  inner  radius  of  the  liner, 
outer  radius  of  the  liner,  specific  gravity  of  rock,  Poisson's  ratio  of 
rock,  Poisson's  ratio  of  the  liner  material,  elastic  modulus  of  rock, 
and  elastic  modtilus  of  the  liner  material. 
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The  streiees  caused  by  the  overburden  may  reach  the  yield 
stresses.  Thus,  the  computer  program  is  organised  so  that  the  first 
computed  static  stresses  are  checked  against  the  yield  stresses  before 
the  dynamic  load  is  brought  into  play. 

Determination  of  Stresses  Caused  by  Ground  Shock.  The  stresses 
caused  by  ground  shock  are  determined  in  two  steps  by  deriving,  first, 
the  impinging  pressure  pulse  and,  second,  the  time -dependent  stress 
distribution. 

Determination  of  Pressure  Pulse.  The  empirically  based  pro- 
ce<hire  in  Reference  4  is  used  to  calculate  the  peak  radial  particle  ac¬ 
celeration,  peak  radial  particle  velocity,  and  peak  radial  displacement 
of  the  free  field  at  the  location  of  the  structure.  According  to  Hugoniot's 
Equation,  the  free-field  stresses  are  approximately  proportional  to  the 
particle  velocity  in  the  elastic  ranger 

O  =  0  •  c  •  V 

^  a 

where  9  s  free-field  stress 
p  s  density 

c^  s  longitudinal  seismic  velocity 
V  =  particle  velocity 

The  plots  of  velocity-time  dependence  and  pressure-time  dependence 
can  be  assumed  to  be  triangles  as  shown  in  Figure  1. 


Figure  1.  Approximation  of  the  preMure-tIme  dependence  by  a  triangle. 
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The  impinging  pressure  pulse  is  fixed  with  these  approximations: 


O  =  o  •  c  •  V 
max  ^  a  max 


2  •  V 

max 

a 

max 


Z  •  d 

max 

d  ^  V 

max 


where  a  =  peak  free  field  pressure 
max 

^max  ~  particle  velocity 
a  =  peak  particle  acceleration 
dmax  ~  radial  displacement 

Magnitude  and  history  of  the  pressure  pulse  are  dependent  upon  the  en¬ 
ergy  coupling  of  the  detonation,  weapon  yield,  energy-conducting  prop¬ 
erty  of  the  rock,  and  distance  from  the  detonation  point.  The  energy 
coupling  of  fully  contained  bursts,  contact  surface  bursts,  and  true  sur¬ 
face  bursts  is  provided  in  the  computer  program.  Empirical  formulas 
for  the  conversion  of  nominal  weapon  yields  into  effective  weapon  yields 
pertaining  to  these  three  kinds  of  detonations  are  given  in  Reference  4. 
The  procedure  for  determining  the  pressure  pulse  is  described  in 
Appendix  B. 

The  energy- conducting  ability  of  the  rock  mass  is  represented  in 
the  empirical  formulas  for  ground  motion  by  the  seismic  velocity.  In 
layered  rock  systems,  the  seismic  velocity  and  the  conductibility  changes 
at  each  interface.  Since  the  incident  wave  and  all  emergent  waves  are 
reflected  and  refracted  at  each  interface,  mathematically  pursuing  the 
complex  energy  flow  is  difficult  and  time  consuming.  Reference  5  pro¬ 
vides  an  approximation  using  an  effective  seismic  velocity  for  the  rock 
body  between  detonation  point  and  structure.  The  effective  seismic 
velocity  is  defined  as  the  quotient  of  the  distance  Rg  and  the  shortest 
transit  time  from  the  point  of  detonation  to  the  structure  (Figure  2): 

R 

s 

C  =  - - 

P  *mi»v 

where  c^  =  effective  seismic  velocity 
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K  -  distance  between  point  of  detonation  and  structure 

S 

t  =  shortest  transit  time, 
min 

A  method  for  determining  the  effective  seismic  velocity  is  described 
in  Appendix  C. 


FIgur*  2.  Layout  for  computing  tho  offoctivo  Miimic  volocity  for  o  loyorod  lyitom. 
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Determination  of  Dynamiolly  Couied  Streai  Diatribution,  Mow 
describes  a  solution  for  computing  the  stress  distribution  caused  by  a 
traveling  harmonic  wave  in  a  liner  of  arbitrary  thickness  and  in  the 
medium  around  the  liner.  In  spite  of  the  analysis  being  limited  to  ho> 
mogeneous  and  isotropic  media  strained  within  the  elastic  range.  Mow's 
method  was  found  to  be  accurate  and  promising  for  further  development. 
Therefore,  it  was  used  for  the  computer  program. 

The  previously  mentioned  limitation  to  the  elastic  range  applies 
only  to  the  medium  in  the  vicinity  of  the  opening  (~3  diameters)  and  to 
the  liner  material.  The  computational  treatment  of  the  pulse  propaga¬ 
tion  from  the  point  of  detonation  to  the  vicinity  of  the  structure  is  not 
based  on  the  elastic  behavior  of  the  rock  media. 

Combining  the  choice  of  Mow's  method  with  the  assumption  of 
homogeneous,  isotropic,  and  elastic  media,  less  energy-absorbing  be¬ 
havior  is  indicated  than  in  real  materials.  This  means  that  the  appli¬ 
cation  of  Mow's  solution  regarding  stress  values  is  conservative. 

The  incident  wave  (Figure  3)  is  described  by  its  potentials: 

longitudinal  wave  ^ 
shear  wave  $  =  0 


where  0 

o 

=  amplitude 

U 

=  circular  frequency 

a 

-  u  /c^  =  wave  number 

6 


1(q(x  -ctft) 


X 


Figur*  3.  Tuniwl  lin«r  wUh  Impinging  itrnu  wov*. 

The  incident  weve  is  reflected  and  refracted  at  the  interface  between 
rock  and  liner  and  is  reflected  at  the  inside  of  the  liner.  Each  reflec> 
tion  and  refraction  causes  a  longitudinal  and  a  shear  wave  that  again 
can  be  reflected  and  refracted.  All  these  waves  contribute  to  the 
particular  motions  and  to  the  stress  distribution. 

The  boundary  conditions  that  guarantee  consideration  of  these 
waves  are: 


1.  The  radial  stresses  in  the  interface  between  rock  and  liner 
are  equal, 

2.  The  shear  stresses  in  the  interface  between  rock  and  liner  are 
equal, 

3.  The  radial  displacements  in  the  interface  between  rock  and 
liner  are  equal, 

4.  The  circumferential  displacement  in  the  interface  between 
rock  and  liner  are  equal, 

5.  The  radial  stresses  at  the  inner  face  of  the  liner  are  zero, 
and 

6.  The  shear  stresses  at  the  inner  face  of  the  liner  are  zero. 
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These  six  boundary  conditions  are  sufficient  to  solve  the  problem.  'I’he 
derivation  of  the  stress  equations  is  given  in  Appendix  D. 

The  approximately  triangular -shaped  pressure  pulse  caused  by 
the  ground  shock  can  be  considered  the  sum  of  an  infinite  number  of 
harmonic  waves,  as  shown  in  Appendix  E.  The  nonzero  potential  of 
the  pressure  pulse  can  be  formulated  by  using  the  Fourier  Integral: 

^(x,  t)e  ‘Oft 

•(B 

The  time -dependent  potential  0  (x,  t)  and  the  frequency-dependent 
potential  <(>(0))  can  be  computed  for  the  pressure  pulse  by  a  method 
given  in  Reference  6.  Stresses  in  terms  of  the  potential  are  given  in 
Appendix  F.  Subsequent  retransformation  of  the  stresses  from  the 
cc  -space  to  the  t-space  leads  to  no  major  difficulties; 

ff  (x*  t)  =  /*<»  (x,  0) )  e  *du; 

Superposition  of  Static  and  Dynamic  Stresses.  The  calculation  of 
static  and  dynamic  stresses  results  in  radial  normal  stresses,  circum¬ 
ferential  normal  stresses,  and  shear  stresses.  As  is  permissible  in 
the  linear  elastic  case,  the  superposition  is  accomplished  by  summing 
the  pertaining  stresses: 

r  s  r  static  r  dynamic 

•  ^B  itatic  ^  ^6  dynamic 

T  •  T  4'  T 

r9s  rd  static  r9  dynamic 

For  further  use,  an  advantage  is  having  the  principal  stresses  available. 
For  the  static,  dynamic,  and  superimposed  case,  the  principal  stresses 
are  determined  by  the  equation: 


8 


Determination  of  Rock  Strength  and  Liner  Material 

Mohr's  Theory  of  Failure  was  chosen  as  failure  criteria  for  both 
the  rock  and  the  liner  material.  The  Air  Force  Special  Weapons  Center^ 
and  the  Geolo(>ical  Society  of  America^  published  test  results  proving 
that  the  envelopes  of  Mohr's  circles  for  most  rocks  can  be  considered 
straight  lines.  This  simple  relationship  between  the  state  of  stress  and 
the  strength  of  rock,  known  as  the  Mohr-Coulomb-Navier 's  Theory,  was 
accepted  for  the  computer  program.  The  equation  for  the  allowable 
maximum  principal  stress  is  derived  in  Appendix  G.  The  value  of  the 
allowable  maximum  principal  stress  is  influenced  positively  by  the 
minimum  principal  stress,  the  cohesive  strength  of  the  rock,  and  the  ^ 
angle  of  internal  friction.  The  American  Society  for  Testing  Materials^ 
recommends  a  curved  line  as  Mohr's  envelope  for  concrete.  Follow¬ 
ing  this  proposal  and  information  given  by  the  Concrete  Division  of 
the  Waterways  Experiment  Station  in  Vicksburg,  Miss.,  the  allowable 
maximum  principal  stress  is  taken  as: 

=  h  •  H  ^ 

^allowable  L  M 

where  =  minimum  principal  stress 
N  =  constant  -  1.  37 
L  -  f ^  =  tensile  strength  of  concrete 


c 


f  =  compressive  strength  of  concrete 
c 

In  steel  construction,  the  increase  of  strength  with  increasing  confine¬ 
ment  is  not  used  for  designing  construction  members.  Information 
about  Mohr's  envelope  of  steel  is  given  in  Reference  10.  The  published 
test  resiilts  do  not  show  an  angle  of  internal  friction  for  steel.  The  en¬ 
velope  is  parallel  to  the  major  stress  axis;  steel  is  considered  a  "von 
Mises"  material.  For  computational  convenience,  the  equation  of  Mohr's 
envelope  for  steel  was  written  in  the  form  of  the  corresponding 
equation  for  concrete: 
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N 


allowable 


where  N  =  M  =  1 


L  =  0F  =  yield  atress  of  steel 

y 

D«t«nn{natfon  of  Survival  Distonco 

Since  the  strength  of  the  rock  and  liner  material  depends  on  the 
state  of  stress  in  these  media,  a  strength  distribution  as  well  as  a 
stress  distribution  will  occur  around  the  underground  opening.  There 
fore,  it  is  necessary  to  check  the  state  of  stress  against  the  state  of 
strength  around  the  opening.  Should  stress  exceed  strength  anywhere 
at  anytime,  the  underground  protective  structure  is  considered  a 
violation  of  the  prerequisite  condition  of  elastic  behavior  of  all 
materials.  As  previously  mentioned,  the  check  of  stress  against 
strength  is  made  twice  for  each  point  and  time: 

1.  The  state  of  stress  caused  by  the  overburden  is  checked 
against  the  strength,  and 

2.  The  state  of  stress  caused  by  overburden  and  ground  shock 
is  checked  against  the  strength. 


If  the  first  check  is  negative,  the  calculation  will  be  stopped  before  the 
dynamic  stresses  are  considered.  The  calculation  also  will  be  stopped 
if  the  second  check  is  negative.  A  negative  resxilt  for  this  second  check 
means  that  the  structure  failed  during  ground  shock.  The  computer  pro¬ 
gram  is  organized  in  such  a  manner  that  the  weapon,  geology,  and  struc¬ 
ture  values,  and  the  depth  of  the  structure  are  kept  constant  for  each 
calculation  while  the  horizontal  distance  from  the  detonation  point  to 
the  structure  increases  for  each  loop  as  long  as  the  structure  does  not 
fail.  The  result  of  one  calculation  cycle  is  the  horizontal  survival 
distance  for  a  preselected  structure  with  given  values  of  cross- 
sectional  size,  liner  material,  liner  thickness,  depth,  geology  con¬ 
dition,  and  threat  attack  procedures. 


Determination  of  Survivability 

The  effectiveness  of  a  protective  structure  is  denoted  by  the 
single -shot  survivability.  Single -shot  survivability  means  the  probability 
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of  survival  of  a  protective  structure  undergoing  a  single  attack. 

For  this  determination  of  survivability,  the  entire  structure  is 
assumed  to  fail  even  if  failure  will  occur  only  in  one  of  its  cross  sec¬ 
tions.  The  survivability  is  influenced  by  weapon  values  and  system 
parameters.  The  main  weapon  values  are  the  hit  accuracy  and  the 
damage  range.  The  hit  accuracy  is  introduced  into  the  computation  of 
the  survivability  by  the  standard  deviations  of  the  weapons  system  in  the 
longitudinal  and  transverse  direction,  respectively.  If  the  standard  de¬ 
viations  are  the  same  in  both  directions,  the  accuracy  of  the  weapon 
system  is  expressed  by  the  Circular  Error  Probable  (CEP).  The  re¬ 
lationship  between  standard  deviation  and  CEP  is: 

9  =  0.  84932  CEP 

The  horizontal  damage  range  of  a  weapon  for  a  particular  structure 
equals  the  horizontal  survival  distance  of  the  structure  calculated  in 
the  design  computation.  The  horizontal  survival  distance  R  determines 
the  aim  area  of  the  underground  structure  at  the  surface  (Figure  4). 

•action  plan 


If  the  design  weapon  hits  inside  this  area,  the  structure  will  fail.  The 
calculation  of  the  survivability  is  derived  in  Appendix  H.  For  this  der¬ 
ivation,  References  11,  12,  and  13  were  used. 

Determination  of  the  Grom  Construction  Costs 
ond  Costs  Per  Usable  Volume  Unit 

Information  about  estimating  costs  for  tunnels  in  rock  was  given 
by  the  Department  of  Water  Resources  of  the  State  of  California^^  and 
the  Metropolitan  Water  District  of  Southern  California  (unpublished 
guidelines).  The  information  given  by  the  Department  of  Water  Resources 
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is  based  upon  data  compiled  from  99  tunnel  projects  around  the  world. 
For  the  most  part,  this  information  was  used  in  the  computer  program. 
The  guidelines  of  the  Metropolitan  Water  District  were  used  as  basic 
information  about  the  costs  of  steel  liner  and  for  checking  and  updating 
the  figures  furnished  by  the  Department  of  Water  Resources. 

The  cost  estimation  procedure  is  due  specifically  to  conventional 
tunneling  methods.  Under  certain  limitations,  machine  excavating 
costs  for  tunneling  are  estimated  at  60%  of  the  costs  of  conventional  ex¬ 
cavating.  This  estimate  is  based  on  information  given  in  Reference  15. 

The  total  construction  costs  are  the  sum  of  the  costs  for  a  main 
tunnel,  access  tunnels,  and  access  shafts.  The  total  cost  for  each 
structure  comprises  liner  costs  and  excavation  costs.  The  liner  costs 
are  calculated  by  "unit  costs  in  place.  "  The  excavation  costs  are  spec¬ 
ified  in  costs  of  labor,  equipment,  energy,  explosives,  drill  bits,  and 
rods.  The  specified  costs  are  classified  according  to  the  size  of  tunnel 
sections  and  according  to  rock  strength.  The  rock  strength  has  been 
described  thus: 

Type  1:  Dry,  massive,  moderately  jointed  or  dry  intact  rock 

Type  2:  Dry,  stratified  or  schistose  rock 

Type  3;  Dry,  moderately  blocky  and  seamy  rock 

Type  4:  Dry,  very  blocky  and  seamy  rock. 

This  classification  has  to  be  expressed  numerically  to  be  used  in  the 
computer  program.  The  approximate  numerical  classification  assumed 
for  the  program  is  based  on  the  cohesive  strength  of  rock: 

2 

Type  1:  K  a  200  kp/cm 

Type  2:  200  kp/cm^  >  K  a  100  kp/cm^ 

2  2 

Type  3;  100  kp/cm  >  K  a  50 kp/cm 

2 

Type  4:  50  kp/cm  >  K 

Besides  the  variables  of  tunnel  length  and  radius,  the  provision  of  ma¬ 
chine  tunneling  depends  upon  the  kind  of  rock.  The  limits  ^  ^  of  machine 
tunneling  are  assumed  as: 
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2 

180  kp/cn»  >  K  >60  kp/cn/ 

L  >  ^000  m 
b  <  5  Ill 

where  L  =  len(;th  of  the  tunnel 

b  =  outer  radius  of  the  tunnel 

Comparable  unit  construction  costs  are  obtained  by  dividing;  the  total 
construction  costs  by  the  usable  volume  of  the  main  tunnel.  A  measure 
of  the  cost  effectiveness  of  the  protective  structures  is  the  quotient  of 
unit  construction  costs  per  unit  of  survivability: 

C  _ _ 

eff  V  .  S 
us 

where  C  =  cost  effectiveness 
eff 

C  =  total  construction  costs 

V  =  usable  volume 

us 

S  =  survivability 

Since  n  ~  1/S,  the  number  of  structures  necessary  to  provide  100% 
survivability  for  one  of  them,  the  cost  effectiveness  quotient  can  be 
interpreted  as  the  costs  for  sure  protection  of  one  unit.  The  cost 
computation  procedure  is  derived  in  Appendix  I.  In  all  the  cost 
equations,  an  adjustment  factor  is  provided  for  updating  and  trans¬ 
forming  into  foreign  currencies. 

COMPUTER  PROGRAM 

Many  portions  of  the  analysis  do  not  require  the  effective  seis¬ 
mic  velocity  between  detonation  point  and  structure  location.  So,  in  the 
interest  of  economy  and  clarity,  two  separate  programs  were  written. 
The  main  program  computes  the  survival  distance,  the  survivability, 
and  the  economic  features.  The  second  program  determines  the  ef¬ 
fective  seismic  velocity  between  detonation  point  and  structure  loca¬ 
tion  in  a  layered  rock  system.  The  second  program  can  be  introduced 
easily  into  the  main  program  as  a  subroutine. 
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Both  programs  are  written  in  FORTRAN  IV  language.  The  time 
necessary  for  one  run  of  the  main  program  using  the  CDC  6600  is  40 
seconds.  The  description,  the  flow  chart,  and  the  listing  of  the  main 
program  is  given  in  Appendix  J.  The  flow  chart  and  the  listing  of  the 
program  for  computing  the  effective  seismic  velocity  of  a  layered  rock 
system  is  given  in  Appendix  K. 

DETERMINATION  OF  OPTIMUM  DESIGN  CONDITIONS 
Example  of  Protective  Structure 

The  total  structure  consists  of  the  main  tunnel,  the  two  access 
tunnels,  and  the  two  shafts  (Figure  5)  with  lengths: 


Length  of  the  main  tunnel  100.  0  m 

Total  length  of  the  access  tunnels  100.0  m 

Inner  radius  of  the  access  tunnels  Z.  0  m 

Inner  radius  of  the  shafts  Z.  0  m 

Inner  radius  of  the  main  tunnel  Variable 


I 


1/ 


r 


T" 

main  tunrrel 


access  tunnel 


Figure  5.  Configuration  sample  for  an  underground  protective  structure  in  rock. 


The  main  tunnel,  access  tunnels,  and  access  shafts  are  lined  by  the 
same  material,  either  reinforced  concrete  or  steel.  The  liner  thick¬ 
ness  T^  of  the  access  tunnels  is: 


T 

a 


T  . 
m 


R 

a 


and  at  least  15.0  centimeters  for  concrete  of  1.Z7  centimeters  for 
steel  where  T^  =  liner  thickness  of  the  main  tunnel 
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=  inner  radius  of  the  main  tunnel 
Rjj  =  inner  radius  of  the  access  tunnels 

The  thickness  of  the  shaft  T,  liner  is: 


T 


s 


1.5  •  T 


a 


and  at  least  20.0  centimeters  for  concrete  or  1.91  centimeters  for  steel.  These 
determinations  are  based  on  rules  of  conventional  mining  and  engineering.  For 
reinforced  concrete  liner,  these  properties  are  auumed. 

350  kp/cm^  (5.000  psi) 


Q)mprestive  strength 
Reinforcement  ratio 
Tensile  strength 
Specific  gravity 
Poisson's  ratio 
Elastic  modulus 

The  properties  for  the  steel  liner  are: 
Yield  stress 
Specific  gravity 
Poisson's  ratio 
EUttic  modulus 


1% 

50  kp/cm^ 

2.4 
0.19 

3.5  •  10®  kp/ci 


m' 


2.52  •  10®  kp/cm2  (36,000  psi) 

7.85 

0.25 

2.05  •  10®  kp/cm2 


Two  conditions  for  the  rock  were  assumed: 


1.  The  entire  rock  body  between  detonation  point  and  structure 
location  consists  of  the  same  rock  type. 

2.  The  earth  body  under  consideration  consists  of  different  rock 
types  or  upper  layers  of  soil. 

For  both  cases,  a  particular  sandstone  and  a  particular  granite  were  selected 
u  the  directly  surrounding  structure.  The  properties  of  the  sandstone  are: 

Specific  gravity  2.5 

Poisson’s  ratio  0.15 

Elastic  modulus  3.0  '  10®  kp/cm^ 
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Longitudinal  seismic  velocity 

3,525  msec 

Factor  of  internal  frietbn 

1.5 

Cohesive  strength 

100  kp/cm^ 

Tensile  strength 

30  kp/cm^ 

of  the  granite  are: 

Specific  gravity 

2.6 

Poisson’s  ratio 

0.20 

Elastic  modulus 

4.0  •  10^  kp/cm 

Longitudinal  seismic  velocity 

4.100  m/s 

Factor  of  internal  friction 

1.5 

Cohesive  strength 

140  kp/cm^ 

Tensile  strength 

60  kp/cm^ 

The  detoiution  ia  assumed  to  be  s  true  surface  burst.  Computations  are  made 
for  100-  and  1,000-kiloton  weapon  yields.  The  protective  structure  is  considered  to  be 
the  target  itself;  the  middle  of  its  target  area  is  Desired  Ground  Zero.  The  CEP  of  the 
attacking  system  it  assumed  to  be  402.2S  meters  (quarter  of  a  mile).  In  subsequent 
comparisons,  the  influence  of  rock  type,  liner  material,  liner  thickness,  and  cross-sectional 
rise  are  presented. 


Optimum  Design  Conditiont  in  a  Rock  Body  of  the  Same  Rock  Type 

Protective  Structures  with  Concrete  Liner  in  Sandstone.  For  an  inner  radius  of 
the  main  tunnel  of  6  meters  and  a  weapon  yield  of  100  kilotons,  the  concrete  liner 
thickneu  was  varied  from  0.5  to  3  meters  as  a  function  of  survivability.  The  surviva¬ 
bility,  which  is  dependent  upon  the  structure  depd^,  is  shown  in  Figure  6.  The 
incr^ng  liner  thickness  doa  not  affect  tignificut|y  the  survival  diKance  for  values 
between  110  and  120  meters  and  a  survivability  of  about  92%. 

The  curves  of  Figure  6  show  that  only  in  tfhe  depth  range  from  50  to  100 
meters  does  the  survival  distance  decrease  wi^  an  increase  of  depth.  Further  increases 
in  depths  beyond  100  meters  are  not  advantageous.  Figure  7  illuNrates  the  total  costs 
for  protection.  The  lower  parts  of  the  curves  represent  the  depth  range  from  50  to 
100  meters.  Increased  eapmiditure  in  this  range  results  in  more  protection,  in  a  decrease 
of  survival  distance,  and  in  an  increaK  of  survivability. 


IB 


MD  KT  1000  KT 
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FIgur*  6  (C).  Survival  dlttonc*  in  a  rock  body  of  tho  tamo  rock  typo  (U). 
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Totol  Costs  (10”  $) 


CONFIDENTIAL 


Inner  Radius  of  Structures;  6  m 
Weapon  Yield:  l(X)  KT 


Survival  Distance  (m) 


Survivability  Rote  (%) 


Note:  Survivability  for  CEP  =  402,25  m 

Figure  7  (C).  Total  costs  for  concrete  structures  in  sandstone  1 

with  different  liner  thicknesses  (U).  * 

I 
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Obviously,  the  slope  of  the  curve  for  the  liner  thickness  of  0.  5 
meter  is  less  than  the  others.  This  means; 

1.  For  the  same  amount  of  additional  money,  the  structure  with 
the  smaller  liner  thickness  offers  more  additional  protection. 

2.  Planning  concrete  structures  with  a  liner  thickness  of  0.  5 
meters  for  a  survivability  of  50%  is  not  economically 
reasonable  since  a  relatively  small  increase  in  the  cost  in¬ 
creases  the  survivability  to  91%- 

Figure  8  exhibits  the  cost  per  usable  volume  and  the  costs  per 
usable  volume  related  to  the  survivability.  The  difference  between 
these  two  cost  values  is  smallest  for  the  liner  thickness  of  0.  5  meter 
and  thus  demonstrates  the  economic  superiority  of  this  design  feature 
as  compared  to  the  others.  Considering  the  previously  mentioned  ad¬ 
vantages,  only  the  liner  thickness  of  0.  5  meter  was  chosen  for  a  com¬ 
parative  study  of  the  cross-sectional  sizes  of  concrete  structures 
having  inner  radii  of  3,  6,  and  9  meters. 

Figure  6  shows  that  within  the  range  considered,  the  section 
radius  does  not  have  a  substantial  influence  on  the  survival  distance 
and  survivability.  The  figure  also  shows  that  the  optimum  depth 
(optimum  in  relation  to  the  survival  distance  and  survivability)  changes 
only  slightly  with  the  cross-sectional  size  (the  optimum  depth  values  for 
the  inner  radii  of  6  and  3  meters  can  be  considered  the  same).  The  op¬ 
timum  depth  for  the  inner  radius  of  9  meters  is  about  120  meters  for  a 
weapon  yield  of  100  kilotons  and  about  220  meters  for  a  weapon  yield  of 
slightly  1.000  kilotons  deeper  than  the  others  having  a  depth  of  100  and 
200  meters,  respectively.  Figure  6  makes  numerically  comprehensible 
the  influence  of  ^e  weapon  yield  on  the  survival  distance  and  surviva¬ 
bility.  The  survival  distance  of  concrete  structures  increases  from  120 
meters  for  100  kilotons  to  220  meters  for  1,  000  kilotons.  Assuming 
the  CEP  equals  402.  25  meters,  the  survivability  decreases  from  about 
92  to  77%.  This  means  a  change  from  a  high  to  a  medium  survivability. 
The  difference  in  the  survivability  becomes  even  larger  when  the  CEP 
of  the  attacking  system  is  smaller. 

Figure  6  illustrates  that  the  larger  weapon  yields  require  a 
greater  depth  for  optimum  survival  distance.  The  economic  aspects 
of  varying  the  cross-sectional  size  is  illustrated  in  Figures  9  and  10. 
Considering  the  total  cost,  the  structure  with  the  inner  radius  of  3 


19 


CONFIDENTIAL 


Costs  per  Us^l*  Volume  (S/m  ) 
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Inner  Rodius  of  Structures:  6  m 
Weapon  Yield:  100  KT 


94  92  90  .  85  80 

Survivability  Rote  (%) 

Note:  Survivability  for  CEP  =  402.25  m 

Figure  8  (C).  Costs  per  usable  volume  for  concrete  structures  in 
sorrdstone  with  different  liner  thicknesses  (U). 
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Concr»f«  Liner  of  0.5  m  ihicknen 
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Concrate  Liner  of  0.5  m  thickness 
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120'^160'W200 

Survival  Distance  (tn) 


240**^280 


Survivability  Rote  (%) 

Note:  Survivability  for  CEP  ■  402.25  m 

Figure  10  (C).  Costs  per  usable  volume  for  structures  in  sandstone 
with  different  inner  radii  (U). 
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meters  is  least  expensive;  however,  structures  of  larger  radii  yield 
favorable  costs  per  usable  volume.  These  economic  considerations 
indicate; 

1.  Small  cross-sectional  sizes  should  be  preferred  for  the  pro¬ 
tection  of  installations  with  small  space  requirements. 

2.  For  volume- oriented  installations  (for  instance,  as  depots),  a 
smaller  number  of  structures  of  larger  cross-sectional  size 
is,  even  considering  survivability,  more  economical  than  a 
larger  number  of  structures  of  smaller  cross-sectional  size. 

Protective  Structure  with  Steel  Liner  in  Sandstone.  The  steel 
liner  is  assumed  to  be  a  tube  of  equal  thickness.  And  the  space  between 
the  rock  and  steel  liner  is  filled  with  concrete.  For  the  scope  of  this 
study,  consideration  of  a  liner  thickness  of  2.  54  cdntimeters  (1  inch) 
was  sufficient. 

Figure  6  shows  the  survival  distance  is  dependent  upon  the  struc¬ 
ture  depth  and  the  weapon  yield.  For  both  the  cross-sectional  sizes 
with  inner  radii  of  9  and  6  meters,  the  survival  distance  decreases 
with  increasing  depth  to  the  optimum  depth  of  150  meters  for  a  weapon 
yield  of  100  kilotons  and  to  250  meters  for  a  weapon  yield  of  1,  000 
kilotons.  The  cross-sectional  size  with  an  inner  radius  of  3  meters 
has  an  optimum  depth  of  about  250  meters  for  100  kilotons  and  offers 
a  smaller  survival  distance  and  higher  survivability.  This  size  effect, 
however,  is  lost  for  the  1,  000-kiloton  weapon  yield.  With  the  higher 
values,  the  weapoh  yield  becomes  the  more  dominant  influencing  factor. 

The  comparison  between  concrete  and  steel  liner  demonstrates 
the  superiority  of  steel  liner.  The  survival  distance  for  the  considered 
steel  liner  is  about  50  meters  less  than  that  for  a  concrete  liner.  The 
superiority  of  a  steel  liner  is  corroborated  by  its  economical  features, 
as  given  in  Figures  11  and  12.  In  spite  of  the  higher  material  unit  costs, 
steel  liner  is  less  expensive  than  concrete  (Figures  13,  14,  and  15)  be¬ 
cause  expensive  excavation  costs  can  be  reduced. 

Protective  Structure  with  Steel  Liner  in  Granite.  The  curves  for 
the  survival  distance  of  structures  in  granite  with  a  steel  liner  (Figure 
6)  show  that  full  protection  against  a  100-kiloton  weapon  yield  can  be 
provided  for  the  three  considered  cross-sectional  sizes  at  a  depth  of 
about  190  meters.  The  surrounding  granite  provides  better  protection 
than  sandstone  from  a  100-kiloton  weapon.  However,  the  superiority 
of  granite  to  sandstone  almost  vanishes  during  a  shock  wave  of  a  1,  000- 
kiloton  weapon.  The  optimum  survival  distance  is  about  160  meters. 
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SImI  Llfwr  of  2.54  cm  thick 
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Figure  11  (C).  Total  costs  for  structures  in  sondsteno  with  differont  inner  radii  (U). 


Volume  (S/m  ) 


CONFIDENTIAL 


Note:  Survivability  for  CEP  =  402.25  m 

Figure  12  (C).  Costs  per  usable  volume  for  structures  in  sandstone 
with  different  inner  radii  (U). 
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Concrvt*  of  0.5  m  thickrms 
Stool  of  2.54  cm  thicknou 
inner  rodiut  of  ftructuro*:  9.0  m 
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Survivobility  Rate  (%) 

Note:  Survivability  for  CEP  =  402.25  m 

Figure  14  (C).  Total  costs  for  structures  m  sandstone  with  steel  and  concrete  liner  (U). 
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The  optimum  depth  fur  granite  is  slightly  less  than  that  for  sandstone 
The  economic  features  of  structures  in  granite  with  a  steel  liner  are 
shown  in  Figures  16  and  17.  They  emphasise  that  full  protection  can 
be  obtained  against  a  100-kiloton  weapon  yield  with  the  same  amount 
of  money  that  is  necessary  to  provide  high,  but  limited,  protection  in 
sandstone. 


Optimum  Design  Conditions  in  o  Rock  Body 
Consisting  of  Different  Types  of  Rocks 

Structures  are  assumed  in  this  part  of  the  study  to  be  embedded 
in  either  the  same  sandstone  or  granite  as  in  the  previous  section.  How¬ 
ever,  the  rock  faces  make  changes  in  the  rock  body  between  the  detona¬ 
tion  point  and  the  structure  location.  Figure  18  illustrates  two  possibili¬ 
ties  of  such  structures. 

Two  cases  are  investigated:  a  structure  of  an  inner  radius  of  6 
meters  embedded  in:  (1)  sandstone  with  a  concrete  liner  and  (2) 
granite  with  a  steel  liner.  For  both  cases,  the  effective  longitudinal 
seismic  velocity  between  detonation  point  and  structure  location  is 
varied  in  1 , 000-meter-per-second  steps  from  1,000  to  6,000  meters 
oer  second. 

Protective  Structure  with  Concrete  Liner  in  Sandstone.  The  sur¬ 
vival  distance  and  the  survivability  are  graphed  in  Figure  19.  For  the 
weapon  values  of  100  and  1,000  kilotons,  the  survival  distance  increases 
with  increasingly  effective  seismic  velocity.  Accordingly,  the  survivabil¬ 
ity  decreases  with  increasingly  effective  seismic  velocity.  Full  protec¬ 
tion  is  given  for  an  effective  seismic  velocity  of  1,000  meters  per  sec¬ 
ond  and  a  weapon  yield  of  100  kilotons.  The  optimum  structure  depth 
increases  slightly  with  dec reasingly  effective  seismic  velocity.  This 
increase  of  the  optimum  depth  is  more  distinct  in  Figure  6  for  the  case 
with  a  1 , 000-kiloton  weapon  yield. 

Figures  20,  21,  and  22  show  the  economic  features.  The  curves 
demonstrate  that  the  expenditures  for  protective  structures  in  a  rock 
body  of  different  types  of  rocks  produce  greater  protection  if  the  effec¬ 
tive  seismic  velocity  is  small.  This  result  shows  numerically  the  sig¬ 
nificance  of  proper  site  selection. 

In  the  case  of  a  1, 000-kiloton  weapon  yield,  the  curves  of  the 
costs  per  usable  volume  as  related  to  the  survivability  have  minimum 
values  for  the  seismic  velocities  of  2,000,  4,000,  5,000,  and  6,000 
meters  per  second.  The  cost  effectiveness  is  the  highest  for  the  sur¬ 
vival  distances  and  depths  pertaining  to  these  minimum  values. 
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Steel  Liner  of  2.54  on  thickness 
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Survivability  Rate  (%) 

Note:  Survivability  for  CEP  -  402.25  m 

Figure  16  (C).  Total  costs  for  structures  in  granite  with  different  inner  rodii  (U). 


Costs  por  Usable  Volume  (S/m^ 
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Survival  Distance  (m) 

L  J  ..j. _ I _ .1  I.  1 - 1 - 1 - 1 

100  9»98  «  M  «90  M  80  75 

Survivability  Rote  (%^ 

Note:  Survivability  for  CEP  *  402.25  m 

Figure  17  (C).  Costs  per  usable  volume  for  structures  in  granite 
with  different  inner  radii  (U). 
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FIgurt  18.  Roek  bodtM  eomltting  of  dtfforont  typos  of  rock. 


Protective  Structure  with  Steel  Liner  to  Granite.  Figure  1 9 
•howe  that  full  protection  can  be  reached  againet  a  100>kiloton  weapon 
yield  for  all  values  of  the  effective  seismic  velocity.  However,  a  small 
seismic  velocity  is  advantageous  here.  too.  because  full  protection  can 
be  realized  for  less  depth  and.  therefore,  for  less  cost.  In  the  case  of 
a  1. 000 -kiloton  weapon  yield  (Figure  19).  full  protection  is  given  for 
values  of  the  effective  seismic  velocity  of  1.  000  and  2.  000  meters  per 
second.  The  curves  for  the  other  values  of  the  effective  seismic  vel¬ 
ocity  show  more  explicitly  that  the  optimum  depth  decreases  with  in¬ 
creasingly  effective  seismic  velocity.  The  cost  curves  in  Figures  20, 
21,  and  22  show  that  structures  in  granite  with  steel  liners  are  more 
economical  than  structures  in  sandstone  with  concrete  liners. 

FINDINGS 

1.  Deep  underground  protective  structures  can  be  designed  for 
high  degrees  of  survivability  and  even  for  full  protection. 

2.  Steel  is  a  better  liner  material  than  concrete  for  protective 
structures  in  rock.  A  liner  2.  54  centimeters  (1  inch)  thick 
satisfies  the  requirements  for  protection  against  a  true  sur¬ 
face  burst  of  1. 000  kilotons. 
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Sandstone  Concrete  Liner  Inner  Radius:  6.0  m 
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Figure  19  (C).  Survival  distance  in  o  rock  body  with  different  rock  types  (U). 


Inner  Rodius  of  Structures:  6.0  m 
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Note:  Survivf^ility  for  CEP  =  402.25  m 

Figure  21  (C).  Total  costs  for  structures  in  sandstone  and  granite  in  coses  of  different 
rock  types  between  detonotion  point  and  surrounding  structure  (U). 
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Figura  22  (C).  Cosh  por  usoblo  volumt  for  structuros  in  londston*  and  gronitt  in 
difforont  rock  typos  botwoon  dotonation  point  and  surrounding  stri 
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(^lu/i)  MMIOA  •I9»»n  J»d  I4»03 


1  londtton*  and  gfonlta  In  cosm  of 
I  point  and  surrounding  strueturo  (U). 
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3.  The  weapon  yield  is  the  dotninaiil  factor  influeiicin(j;  the  survi¬ 
val  distance  and  optinuuii  depth.  Its  influence  was  flete r mined 
numerically  for  the  100-  and  1 , 000-kiloton  values.  The  results 
are  shown  in  Figures  6  and  19. 

4.  The  influence  of  the  type  of  enveloping  rock  appears  to  he  more 
significant  for  the  smaller  weapon  yield  of  100  kilotons. 

5.  The  rock  interfaces  between  the  detonation  point  and  the  struc¬ 
ture  are  of  great  importance  since  the  survivability  increases 
greatly  when  the  effective  seismic  velocity  of  the  rock  body  de¬ 
creases.  The  influence  of  the  rock  interfaces,  represented  here 
by  the  seismic  velocity,  is  shown  nunierically  in  Figure  19. 

6.  No  significant  influence  of  the  cross-sectional  size  on  the  survi¬ 
vability  could  be  found  for  structures  with  concrete  liners. 

7.  For  structures  with  concrete  liners,  the  increase  in  resistance 
of  liner  thicknesses  greater  than  0.5  meter  was  found  to  be 
insignificant. 

8.  With  a  steel  liner,  tunnels  having  an  inner  radius  of  3  meters 
provided  greater  survivability  than  those  with  inner  radii  of 

6  and  9  meters. 

9.  Steel  liners  provide  greater  survivability  than  concrete  liners. 

10.  Structures  with  steel  liners  are  less  expensive  than  those  with 
concrete  liners. 

CONCLUSIONS 

From  the  foregoing  study  of  deep  underground  tunnels  in  rock,  it 
can  be  concluded  that: 

1.  An  increase  in  the  depth  of  cover  does  not  always  lead  to  an 
increase  in  survivability.  There  are  optimum  depths  which, 
if  exceeded,  will  increase  the  survival  distance  and  decrease 
the  survivability. 

2.  The  optimum  site  profile  consists  of  a  thick  soil  layer  over  a 
strong  basement  rock. 

The  computer  program  described  herein  gives  the  location,  with  respect 
to  ground  zero,  of  a  tunnel  with  a  preselected  survivability.  It  also  pro¬ 
vides  unit  and  total  costs. 
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APPENDIX  A 


STRESS  DISTRIBUTION  CAUSED  BY  THE  OVERBURDEN 


Ftgur*  23.  Tunnel  eonflguratton. 


Nototiont 

H  s  depth  of  the  structure  in  meters 
m  »  inner  radius  of  the  liner  in  meters 

b  =  outer  radius  of  the  liner  in  meters 

r  s  coordinate  radius  in  meters 
6  =  coordinate  angle  in  degrees 

3 

y  •  specific  weight  of  the  overburden  in  Mp/m 

2 

Ej  ~  elastic  modulus  of  surrounding  rock  in  kp/m 

2 

E2  =  elastic  modulus  of  liner  material  in  kp/m 
Uj  X  Pole  son's  ratio  of  surrounding  rock 
1^2  ■  Poisson's  ratio  of  liner  material 
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®r(2)tt 

’'re(2)«t 


2 

radial  atresa  in  rock  cauaed  by  the  overburden  in  kp/m 

circumferential  atreaa  in  rock  cauaed  by  the  overburden 
in  kp/m^ 

2 

ahear  atreaa  in  rock  cauaed  by  the  overburden  in  kp/m 

2 

radial  atreaa  in  liner  cauaed  by  the  overburden  in  kp/m 

circumferential  atreaa  in  liner  cauaed  by  the  overburden 
in  kp/m^ 

2 

ahear  atreaa  in  liner  cauaed  by  the  overburden  in  kp/m 


The  apecific  weight  ia  the  average  apecific  weight  divided  by  the 
entire  depth  H.  The  elaatic  propertiea  E|  and  of  the  aurrounding 
rock  are  the  reapective  valuea  of  the  rock  body,  within  3  cavity  diam- 
etera,  around  the  atructure.  The  atreaa  equationa  are: 
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APPENDIX  B 


DETERMINATION  OF  THE  PRESSURE  PULSE 


Figur*  24.  Tunnal  location  with  ratpoct  to  grotjnd  zoro. 


Notation: 

H  =  depth  of  the  structure  in  meters 

R  =  horizontal  distance  between  detonation  point  and  structure  in 
meters 

R  =  slant  distance  between  detonation  point  and  structure  in 
meters 

C  =  effective  longitudinal  seismic  velocity  for  the  rock  body 

^  between  the  point  of  detonation  and  the  structure  in  meters 

per  second 

C  =  c  =  longitudinal  seismic  velocity  of  the  rock  body  of  three  tunnel 
^  ^  diameters  around  the  structure  in  meters  per  second 


^max  ~  rock  field  peak  radial  particle  acceleration  in  g's 


V 

max 


free  rock  field  peak  radial  particle  velocity  in  meters  per 
second 


'^max  ~  rock  field  peak  radial  displacement  in  meters 


45 


=  ris*  time  of  the  preeeure  pulse  in  seconds 

t .  =  duration  of  the  pressure  pulse  in  seconds 

d 

t.  =  decay  time  of  the  pressure  pulse  in  seconds 
D 

=  nominal  weapon  yield  in  kilotons 

Three  kinds  of  atomic  weapon  bursts  are  distinguished  and  con¬ 
sidered.  The  fully  contained  burst  is  assumed  to  have  an  energy 
coupling  for  the  ground  shock  of  100%.  The  effective  weapon  yield 
We  ie  equal  to  the  nominal  weapon  yield  W|^: 

W  =  W 

e  n 

The  contact  surface  burst  detonates  a  small  but  finite  distance  above 
the  ground  surface: 

W  =  W  (0.02  -  ) 

e  n  S  n 

The  true  surface  burst  detonates  coincident  with  the  ground  surface: 

W  s  0.01  W  (6.00  -  log  W  )  . 
e  n  n 


The  true  distance  between  the  detonation  point  and  the  center  of  the 
structure  is: 

The  true  distance  is  converted  into  the  scaled  distance: 

R 

The  scaled  distance  is  the  variable  in  the  equations  for  free- 
field  peak  radial  particle  acceleration,  velocity,  and  displacement. 
These  equations  are  obtained  from  the  graphs  in  Reference  4: 
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where  an^axj  dep®"**®  ^e* 


c 

R 

log  a 
^  max. 

1 

p 

e 

c  < 

1800  m/s 

R  a  liO  m 

-4.4064  log  R^  +  9.  1425 

p 

c  « 

1800  m/s 

e 

R  >  130  m 

-1.9991  log  R^  +  4.  0456 

p 

1 800  m/a  <  C  ^ 

D 

3600  m/s 

R  <  200  m 
e 

-3.  9798  log  R  +  9.3843 

r 

1800  m/s  <  C  « 
P 

3600  m/ s 

R  >  200  m 
e 

-1 , 8143  log  R  +  4. 4104 
e 

3600  m/s  <  C 

D 

R  s  105  m 
e 

-2.6394  log  R  +  6. 9549 
e 

tr 

3600  m/s  <  C 

R  >  105  m 
e 

-4. 8804  log  R  +  11. 4999 
e 

P 
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I 


Peak  radUl  particle  velocity. 


V 


V 

max 


max, 

[r 


2440 

V 

max^ 

4880 


where  V 

max 


C 

-E 


c  < 

p 

c  < 

p 

1800  m/e  <  e 
P 

1800  m/e  <  < 

P 

3600  m/s  <  Cp 
3600  m/e  <  C 

P 


(c  ) 

p 

if 

C  s  1800  m/s 

P 

■<s> 

if  1800  m/e 

<  C  a  3600  m/e 

P 

■s> 

if  3600  m/e 

<  C  . 

P 

depende  upon  C  and  R  , 
1  " 

thue: 

p 

log  V 
®  max^ 

l\ 

c 

1800  m/e 

R  £  105  m 

e 

.3.6860  log  R  +  6.8244 

1800  m/e 

R  >  105  m 
e 

-1.3497  log  R^  +  2.0926 

3600  m/e 

R  s  75  m 

e 

.2.9822  log  R  +  6.0408 

3600  m/e 

R  >  75  m 

e 

.1.6660  log  R  +  3.5528 

R  a  120  m 

e 

.2.4134  log  R  +  5.2966 
c 

R  >  120  m 

a 

.1.2544  log  R^  +  2.8790 
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II 

0 

II 

II 

[] 

a 

a 

0 

ii 

i 

i 

11 

0 

0 

LI 


Th«  £r««-{ield  stresses  are  approximately  proportional  to  the 
particle  velocity  (Hugoniot  Equation)} 

9  s  A.  C  *  V  ,  and 

a 

O  a  p‘C  ‘V 
max  a  niax 

where  o  a  free-field  radial  stress  in  kp/m^  , 
p  a  mass  density  =  •  1000 

For  the  triangular- shaped  load  impulse  (Figure  25),  the  rise 
time  is 


2*v 
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and  the  pulse  duration  is 

2-d 

max 

d  ■  V 

max 

For  further  computations,  it  is  convenient  to  have  the  decay  time 
expressed  as: 


Figure  25.  Streu  and  velocity  pulie. 
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APPENDIX  C 


DETERMINATION  OF  THE  EFFEaiVE  SEISMIC  VELOCITY 
OF  A  LAYERED  ROCK  SYSTEM 


Figur*  26.  DIrtct  and  indirtcf  tayi  from  tho  point 
of  dotonotion  to  tho  tlWcturo. 


Notation: 


H  =  depth  of  the  structure  in  meters 

R  =  horizontal  distance  between  detonation  point  and  structure 

in  meters 

R  =  slant  distance  between  detonation  point  and  structure  in 

meters 

T  =  thickness  of  layer  k  in  meters 

Tjj  =  thickness  of  layer  package  between  point  of  detonation  and 

structure  in  meters 

L  =  maximum  thickness  of  the  layer  packages  penetrated  by  the 

max  •  j-  •*. 

indirect  rays  in  meters 

a  =  acute  angle  between  ray  and  the  normal  of  the  refracting 

plane  in  degrees 

f  =  angle  of  layer  inclination 

6  =  angle  between  horizontal  distance  (R)  and  slant  distance 

(R  )  in  degrees 
s 

=  seismic  velocity  of  layer  k  in  meters  per  second 

C  =  effective  seismic  velocity  for  the  distance  between  point  of 

^  detonation  and  structure  location  in  meters  per  second 

k  =  ordinal  number  of  the  layers 

The  effective  longitudinal  seismic  velocity  is  the  quotient  of 
the  slant  distance  from  the  point  of  detonation  to  the  structure  and  the 
shortest  time  necessary  to  travel  this  distance: 


The  structure  can  be  hit  directly  by  a  pressure  ray  (ray  group  I  through 
layers  6,  7,  8,  and  9)  going  through  the  layers  between  the  point  of 
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i.lctonatiDii  and  the  structure,  indirectly  by  ray  j^roup  11  which  penetrates 
into  layers  above  the  structure  and  returns  after  havinj^  traveled  parallel 
in  one  of  the  layers  above,  and  iiulirectly  by  ray  group  111  which  pene¬ 
trates  into  layers  below  the  detonation  point  and  propagates  to  the 
structure  after  having  traveled  parallel  in  one  of  the  layers  below. 

The  travel  times  of  tlie  direct  ray  and  all  the  indirect  rays  are  com¬ 
puted  and  compared.  By  Snell's  Law,  the  direct  ray  has  to  yield  these 
conditions : 


and 


sin  Ot, 
_ 1^ 

sin  a,  , 
kf  1 


E  Xk 

m 

where 

X 


=  X 


=  R*  cos  €  +  H*  sin  f 
=  '8“k 


The  path  of  the  direct  ray  must  be  found  by  iteration.  The  loop 
is  initiated  by  the  ray  that  starts  into  the  straight-line  direction  from 
the  detonation  point  to  the  structure: 

aj  =  90  -  (e  -  6)  . 


A  direct  pressure  ray  does  not  exist  if  its  path  leaves  the  surface  of 
the  half  space  (Figure  27). 
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The  travel  time  is: 


t 


k=n 


E 

k=iu 


T 


k 


cos  a,  • 
k 


C 

k 


In  the  example  of  Figure  27,  m  =  9  and  n  -  6, 

The  computation  of  the  refraction  angles  ,  since  the  indirect 
rays  start  with  the  one  that  is  90  degrees,  is: 

a  =  90° 

r 

sin  a 

r 

sin  a, 
k 

The  path  length  traveled  parallel  to  the  layers  is: 


=  1 


C 

r 


X 

s 


k=j 

k-m 


k=n 

-  E 

k=  j 


For  ray  group  III2  in  Figure  26,  m  =  10,  j  =  11,  and  n  =  6.  If  an 
assumed  ray's  length,  Xg,  is  negative,  that  ray  never  can  be  the  fastest 
one  and  is  to  be  ignored. 

The  travel  time  for  an  indirect  ray  is: 


k=j 


E 

k=m 


T 


cos  (V 


k_ 

k 


k=n 


+  E 

k=  j 


cos  a. 


C 


k 


For  computational  convenience,  the  detonation  point  and  the  structure 
locations  are  assumed  to  divide  their  particular  layers  by  an  imaginary 
cut  into  two  layers  of  the  same  seismic  velocity. 

To  simplify  the  computation  and  the  use  of  the  program,  the  half 
space  is  considered  to  consist  of  four  areas  as  shown  in  Figure  2  8.  The 
areas  are  limited  by  the  relations; 

Area  A;  €  is  positive  and  0  s  6  <  (90  -  e) 

Area  B:  f  is  positive  and  (90  -  f)  ^  6  <  90 

Area  C;  e  is  negative  and  |  c]  ^  6  ^  90 

Area  D:  f  is  negative  and  0  s  6  <  |e|  . 
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Figure  28.  PartlHon  of  the  half  space  In  areas  related  to 
layers'  Inclination  and  structure  location. 


This  distinction  is  based  on  the  following  sign  convention: 

Distance  R  is  always  positive. 

Depth  H  is  always  positive. 

Location  angle  6  is  always  positive. 

Layer  inclination  c  is  positive  if  the  detonation  point  is  located 
"above  the  layers"  and  is  negative  if  the  detonation  point  is 
located  "below  the  layers." 

"Above  the  layers"  means  that  an  observer  moving  into  the  direction 
of  R  meets  lower,  in  general  earlier  sedimented,  layers.  "Below 
the  layers"  means  that  the  observer  meets  higher,  in  general  later 
sedimented,  layers  (Figure  28). 

With  this  area  and  sign  convention,  the  previously  given  expres¬ 
sions  for  this  distance  X  and  for  the  first  approximate  angle  change  to: 

Area  A:  X  =  R  .•  cos  e  -  H  .  •  sin  c 

A  A 

=  90  -  f 

Area  B:  X  =  R  *  cos  c  -  '  siu  e 

B 

a ,  =  +  c  -  90 

1  o 
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Area  C; 


Area  D: 


X  =  R  •  cos  €  +  H  sin 

a,=  90  f  c  -  6^ 


X  =  R 


D 
=  90 


COB  €  -I  sin 

■  C  + 


€ 


C 


All  the  rays  that  are  supposed  to  hit  the  structure  travel  between  the 
perpendiculars  from  the  point  of  detonation  or  the  location  of  the 
structure  to  the  layers.  A  ray  traveling  beyond  a  perpendicular  cannot 
cross  the  perpendicular  because  the  signs  of  both  sides  of  the  basic 
law: 


sinai^l 


sin  otj 


are  the  same  and,  in  the  case  of  refraction,  the  angles  are  allowed  to 
change  only  from  0  to  90  degrees. 

Because  of  these  relations,  ray  group  II  does  not  exist  for  area 
A.  To  penetrate  the  first  layer  above  the  detonation  point,  the  rays 
would  have  to  travel  beyond  the  perpendicular  in  area  D  and  from  there 
they  cannot  return  to  area  A  by  refraction. 

Besides  the  layer  package  between  the  detonation  point  and  the 
location  of  the  structure,  layers  above  and  below  this  package  are 
significant  in  selecting  the  shortest  travel  time.  The  program  user 
will  want  to  know  how  many  layers  above  and  below  to  consider  (Figure 
29).  The  maximum  thickness  of  the  layer  package  above  may  be 
determined  by  the  assumption  that  the  layers  penetrated  do  have  the 
minimum  seismic  velocity  for  rock  of  1,000  meters  per  second  and 
that  the  layer  in  which  the  ray  travels  parallel  does  have  the  maximum 
seismic  velocity  in  rock  of  6,500  meters  per  second. 

From  Figure  29  it  can  be  seen  that: 

direct  1000 


indirect  1000*  cos  6500  1000*  cos 


sin  aj 


£l  ^  1000 

6500 


0.  1539 


cos  ttj  =  0.  9881 
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Figure  29.  AAaxImum  thickness  of  the  layer  package  above  and  below 
the  directly  penetrated  layers. 


tan  ttj  =  0.  1557 

X  =  R*  cos  c  +  H  sin  e 

X  =  (R'  cos  e  +  H  sin  e)  -L-tantt  -  (L  +  T  )•  tan  Or 

S  1  U  i 

=  (R'  cose  +  Hsine)-  0.1557(2*L  + 

t  =  t 

direct  indirect 
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L  I  T 


Since 


If  0  s 


xTz  i 
>  R  t  M 


D 


1000  1000-0.9881  6500  1000-0.9881 


-t  11^  R-  cos  c  +  H  •  sin  f  0.  987  •  T 


D 


1. 974 


1, 974  -  6. 5 


■  R  •  sin  C  4  H  '  cos  e  =  0 


-R-  sin  e  4  H  •  cos  C  _  ^ 


1. 974  -  6.  5 


H 

tan  c  =  ^ 


1. 974 


L" 

e  s  90 


-R-  cos  C  -  H  sin  C 
1.974  -  6.5 


<  0 


Lg  (f  =  arctan  ^ 


COS  f 


sin  C 


1 


V 


tan  f  4  1 
tan  t _ 


V.  ^ 


tan  f  4  1 


R 


Vr^  + 

H 


max 


R^  4  T 


D 


1.  974  1.974  •  6 


.5  Vr^  4  2 


L  =  0.585  -  R^  4  0.5  •  T 

max  ^ 


L  0.6  •  R„  4  0.5  •  T 

max  o  ^ 
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Parallels  to  the  layers  in  the  distance  Lmax  from  either  the 
detonation  point  or  the  location  of  the  structure  determine  the  lower 
and  upper  boundary  of  the  rock  body  which  has  to  be  taken  into  con¬ 
sideration  (A,  B,  C,  and  D  in  Figure  29).  In  the  demonstration  case 
of  Figure  29,  however,  the  distance  Lmax  on  the  perpendicular  into 
the  upper  layers  exceeds  the  ground  surface.  Therefore,  the  real 
upper  boundary  is  the  ground  surface  itself.  So  layer  T2  is  the 
upper  boundary  layer  and  layer  T  |  0  is  the  lower  boundary  layer. 

The  direct  rays  in  areas  A  and  D  and  the  indirect  rays  of 
group  11  in  areas  B,  C,  and  D  have  to  be  checked  against  leaving  the 
half  space.  For  the  check,  a  critical  Xa  of  each  layer  is  computed 
(Figure  29). 


Figure  30.  Propagation  areas. 

Rays  in  area  A  are  considered  real  if  the  condition: 


r  X  s  x 


ak 


is  met.  The  previously  mentioned  rays  in  areas  B,  C,  and  D  are 
considered  real  if  the  condition: 


Ex,  a  X 


ak 


is  yielded. 

The  input  to  the  program  is  organized  in  such  a  manner  that 
the  layers  of  the  considered  rock  body  are  denoted  by  a  continuous 
row  of  arable  numbers  starting  with  "1"  .  The  upper  boundary 
layer,  generally  the  geologically  youngest  layer,  is  layer  number  "1", 

Summarization  of  the  rules  for  preparing  the  geology  input: 

1.  Determine  the  sign  of  the  inclination  angle. 
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2.  Determine  the  side  boundaries  of  the  rock  body  to  be  con¬ 
sidered  by  drawing  the  perpendiculars  to  the  layers  through 
the  point  of  detonation  and  the  location  of  the  structure. 

3.  Determine  the  upper  and  lower  boundaries  by  parallels  to 
the  layers  in  the  distance  Lmax  0.  6Rs  +  0.5  Tq  from  the 
point  of  detonation  and  the  location  of  the  structure. 

4.  Check  whether  the  ground  surface  substitutes  the  upper 
parallel  as  upper  boundary. 

5.  Divide  the  layers  in  which  detonation  point  and  structure  are 
located  by  an  imaginary  cut  into  two  layers  each. 

6.  Enumerate  the  layers  of  the  rock  body  to  be  considered,  start¬ 
ing  with  "1"  for  the  upper  boundary  layer  (generally  the 
geologically  youngest  layer). 

7.  Order  the  layers  penetrated  by: 

Ray  Group  I  to  the  input  notation  NLF, 

Ray  Group  II  to  the  input  notation  NLBl,  and 

Ray  Group  III  to  the  input  notation  NLB2. 
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I 


I 


i  i 


* 

I 

•  r 

i 


I  T'* 

,  J 


J 


DERIVATION  OF  THE  STRESS  DISTRIBUTION  CAUSED  BY  AN 
INCIDENT  HARMONIC  WAVE 


X 

Figure  31 .  Stress  wave  intercepting  cylinder. 


Notation: 

r  =  coordinate  radius  in  meters 
a  =  inner  radius  of  the  structure  in  meters 

b  =  outer  radius  of  the  structure  in  meters 


=  longitudinal  seismic  velocity  in  meters  per  second 

Cq  =  transverse  seismic  velocity  in  meters  per  second 
P 

2 

a  =  radial  stress  in  kp/m 
r 

2 

=  circumferential  stress  in  kp/m 

2 

T  „  =  shear  stress  in  kp/m 

r0 
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1  as  a  subscript  denotes  the  rock  medium 


2  as  a  subscript  denotes  the  liner. 

STRESS  TENSOR 

The  differential  equation  governing  the  motion  of  an  elastic,  iso¬ 
tropic,  and  homogeneous  medium  is; 

«  • 

-  2-  - 
(X  +  U)VV  -u+jiv  u=pu 

where 


X  =  Lame's  1.  constant  of  the  medium 

p  =  Lame's  2.  constant  of  the  medium 

p  =  density  of  the  medium 

u  =  displacement  vector 

The  displacement  vector  can  be  expressed  by  the  derivatives  of  a 
scalar  potential  and  a  vector  potential: 

^  =  70  + 


The  potentials  satisfy  the  wave  equation; 


Of 


2t  1  T 
7  ^ 

where 

0  =  longitudinal  displacement  potential 
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.  > 


.  i 


^  =  distortional  displacement  potential 

Q  _  jXJ^lA  -  longitudinal  seismic  velocity 
a  ^  0 

-  distortional  seismic  velocity 


The  stress  tensor  is  related  to  the  displacement  vector  by: 


where 


a  =  X(y-u)i  +  u(vu  +  u7) 


I  =  unit  tensor 


For  the  assumed  plane  strain  problem,  the  vector-displacement 
potential  is  directed  into  the  axis  of  the  cylinder; 

j|)  =  e  j/) 

z 

The  vector  potential  is  the  scalar  product  of  the  unit  vector  e  and  the 
magnitude  ij) .  Substituting  the  potentials  into  the  stress  tensor,  the 
stress  equation  in  cylindrical  coordinates  can  be  written  as; 


XV 


rr  =117^^  +2u/—  +--^ 

ag  XV  ^  ^^2  +r  3r  r  9r?,fl\2  99  / 

re  ^*'\r9r30  ^2  30  +  r  9  r  +  ^2  2  2  ^  j 


The  equations  for  the  displacement  are: 


u 


_  5  (b  ^  9  <!) 

~  3r  r  30 
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,  i  _  Aill 

fl  r  90  9  r 


STRESSES 

The  potentials  of  an  incident,  longitudinal,  harmonic  wave  propa¬ 
gating  in  the  x-  direction  are; 


i(ax  -  tot) 


where 


0  " 

0  =0 


0  =  amplitude 

o 


0) 

O 


-  JiL 
C 

a 


circular  frequency 

wave  number  of  the  longitudinal  wave 


The  nonzero  potential  can  be  written  in  cylindrical  coordinates  as: 


0  =0  e‘^“*  "  =  0  J  €  i"  J  iar)  cos  (n0)e'^‘^* 

o  o  Q  n  n 


where 


n 


n 


=  Bessel  function  of  first  kind  of  order  n 
=  1  for  n  =  0 
=  2  for  n  ^  0 


The  incident  wave  is  reflected  at  the  liner.  The  outward- propa¬ 
gating  longitudinal  and  shear  waves  are  determined  by: 


^(1)  "  2)  («!*•)  cos  "Be 

n=0 
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m 

(1) 


ir 

y  11  *  (fl  I  *^6  ‘• 

ti  n  ^  1 

n  U 


The  waves  propagating  through  the  liner  towards  the  inner  boundary  are 
represented  by: 


'f 


(1) 

(2) 


>  M  H  (tv^r)  cos  n0  e 

n  "  "  ^ 

n-O 


li' 


(1) 

(2) 


CO 


I 

n=0 


N  H 
n 


(2) 

n 


sin  nQ  e 


The  refracted  wave  is  reflected  at  the  inner  face  of  the  liner.  The  re¬ 
flected  waves  propagate  outward: 


(R) 

^(2) 


eo 


1 

n=0 


R  H 
n 


(1) 

n 


cos  n0e 


1/) 


(R) 

(2) 


y  S  H  (/}  r)  cos  nge  ^ 
n=0 


where 


n 


Hankel  function  of  first  kind  of  order  n  (for  inward- 
propagating  waves) 

Hankel  function  of  second  kind  of  order  n  (for  outward- 
propagating  waves) 


a 

8 

O) 


=  -ItL 
C 


Of 


wave  number  of  longitudinal  wave 


wave  number  of  shear  wave 


=  circular  frequency 
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A  .  B  ,  M  ,  N  ,  R  ,  and  S  are  coefficients  that  have  to  be  deter- 
n  n  n  n  n  n 

mined  by  using  the  boundary  conditions: 

At  the  interface  between  rock  and  liner. 


r 

=  b 

^r{l) 

"  "r(2) 

^r0(l) 

%(1) 

“r(2) 

'"ed) 

"  '^0(2) 

the  innerface  of  the  liner  has  to  be  stress-free: 


^r{2) 

^0(2) 


=  a 

=  0 
=  0 


Superposition  of  the  effects  of  the  incident  reflected  and  refracted 
waves  respectively  leads  to  the  stress  equations; 


=  Zu^r 


-2 


JJ 


0  €  i"  +  A 

on  1  nr  n  1  nr  n  1  nr  / 


6(1) 


-2  V  t  t^(i 


+A 

nr  n  1  nr 


B 

n  1  nr 
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where 


T  -la 

red)  ^ 


+  A  + 

nr  n  1  nr 


n  1  nr/' 


sin  nfle 


-iu;t 


"^(2)  =^“2'' 


-2 


^‘^+N  D^‘*+R 

nr  n2  nr  n2  nr  n2  nr 


X^cos  nQe 


'e(2) 


2u_  r 


X  (cos  n0e 


T  =  2u  r“^  y(M  -E'‘'  +  N  ,€'‘'  +  R  ,E-''+S  , 

r0(2)  2  ^4<l\  n2  nr  n2  nr  n2  nr  n2 

X  (sin  n0e'^‘*^*)j 


’  £  [m  -  N  D 

n 2  nr  n2 

"1 


"'‘'+R  B'R' 

nr  n2  nr  n2  nr 


.(i) 


(i).,.  ^.(R).o  e(R)\ 

i  nr  / 


iDlr  =("^+  ”-i^l 


„(R)  /  2  ^  1  .2  2)  „(l)  ,  ,  „(U  ,  , 

:°nr  M"  +”-2»2'^)“n  <“2'’ ■“2"”n.  1 '“2'’ 


(i)  2 

=(n  +n|J_^(ajr|.  n„jrj_^_jta|r) 

2<r  ="''^'’'<"'<»2--'-"“2'«n-’l'“2--' 
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„(R)  ,  2 

E  =  (n 
1  nr 


.f'*' 

I  nr 


2  nr 


p(R) 

1  nr 


(R) 

2  nr 


-v 


(R) 

1  nr 


-(n 


(R) 

2  nr 


•(n 


2  nr 


=  “(n 


.(R) 

1  nr 


=  -(n^ 


e(R) 

2  nr 


=  -  n 


=-(n^ 


The  coefficients  A  , 
n 


+  n)  H  (a  ,  r)  -  n«  rH  ,  (fV  ,  *■) 
n  1  1  n-  1  1 

+  n)  H|^‘^0,2r)  -  (o^r) 

+  n  -  J^(Oir)  +  a, j (ttj r) 

22  1  2  2\  (2)  (2) 
tn-a^'-  +2«2'')«„  <“2’'' 

.n.„^.^i^^rV''>,a2r).a,<’,(»2r. 

+  n)  hJ^'*  (^,r)  +  (/9,r) 

+  n)  Hjf'  (J^r)  +  n?2''«n.’l  <?2’'' 

B  ,  M  ,  N  ,  R  ,  and  S  can  be  determined  from: 
n’  n  n  n  n 
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where 


ra.l  Ic.l  =  |b.  ! 
L  ijJ  I  j  I  I  1  ) 


i.(R) 

5(R) 

2  nb 

a^lb 

2  nb 

2  nb 

p(R) 

V\^nh 

^(R) 

''rnb 

2  nb 

^(i) 

2-nb 

p(R) 

2^nb 

2^nb 

«!*>»!.' ’s'” 

"“"S'" 

bHn  (o^b) 

(1)' 

bH)^  (O^b) 

""n ’<%'■' 

BH;,"(a,b)  nnjf' (O^b) 

ftbHf  (^b)  nHj^'^bl  (l2bH;,‘''(%b) 

0 

o 

,5“’ 

d(R) 

2  na 

2  na 

2  na 

2  na 

o 

o 

(R) 

=(R) 

oE 

2  na 

2'  na 

2  na 

2  na 
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The  expressions 


.xO'’ 

J  na 


or 


denote  the  values  of 

(yl 

J  nr 

for  r  =  a  or  r  =  b.  respectively. 

After  determining  the  coefficients,  the  stress  can  be  calculated 
from  the  stress  equation. 
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APPENDIX  E 


DERIVATION  OF  THE  FOURIER-TRANSFORMED  POTENTIAL 
OF  THE  TRIANGULAR  PRESSURE  PULSE 

THE  FOURIER  TRANSFORM  IN  GENERAL 


The  Fourier  transform  f(a))  of  the  pulse  f(t)  is 

f(aj)  =-^1  f(t)e'‘^‘dt 
-00 


Because  the  pulse  is  considered  to  start  on  times  greater  than  zero, 
the  one-sided  Fourier  transform; 


f(W) 


_1_ 

2ff 


i(0t 

e 


dt 


will  be  used. 

If  the  function  £(t)  is  approximated  by  straight-line  segments  in 
the  increments  of  Atj  ,  the  function  can  be  written  as; 

Af.(t) 

for 


t.  ,  <t  i  t. 
J-1  J 


where 


Af.(t)  =  f.(t)  -  f.  ,(t) 
J  J  J-1 
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The  first  derivative  of  the  approximating  function: 


df.(t)  Af.(t) 

_1_ _ 2_ 

dt  £st. 

J 


is  continuous  inside  the  interval  At.  .  At  the  limit  arguments,  t.  ^ 
and  tj  ,  the  first  derivative  is  discontinuous.  The  function  of  the? 
first  derivative  is  a  step  function. 

The  second  derivative  inside  the  interval  is  zero: 


0 

dt^ 

but  is  undefined  in  this  form  for  the  limiting  arguments  of  the  interval. 
The  second  derivative  for  the  limiting  arguments  can  be  determined  by 
the  Dirac  impulse: 


2 

Af(t.)  6  (t  -  t.) 
dt  J  ^ 


if  t  =  t. ,  where 

Af'(t.)  =  — 
J 


dt 


df.(t) 

-■J... 

dt 


Since  the  Dirac  impulse  is  zero  for  arguments  different  from  zero: 


6(t-t.)=  0 
J 


If  t  ^  tj  , 

the  second  derivative  for  the  interval  and  its  limit  on  the  right  may  be 
written: 

f"(t)  =  Af ' (t.)  6  (t  -  t.) 

J  3 

t.  ,  <  t  5  t. 

J  -  1  J 
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The  Fourier  transform  of  f* (t)  is: 


CM 


iwt 


dt 


f*(a)) 


=  •57  Z  J  Af '(..)» (t  -  y  e‘"‘ d. 


Using  the  definition  of  the  Dirac  function: 

f  .‘“'i 

■0  ^ 

the  transform  of  f*  (t)  is  simplified  to: 

,  n  iojt. 

iUu))  =-^1  Af'(t.)c  J 
The  second  derivative  of  the  harmonic  wave  form  is: 

r(w)  =-«^f(w) 

Therefore,  the  transform  of  f(t)  can  finally  be  obtained  as  follows; 


f(w) 


f(a}) 


ilM 

0,2 


n  Af ' (t.)  e 

=  —  T  - ^ - 

j  =  l  -0)2 


io}t 


where,  according  to  prior  definition: 

f..l(t)  -  f  (t)  f  (t)  -  f  (t) 
^  J+1  J  J  J-1 
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APPLICATION  OF  FOURIER  TRANSFORM 


The  free- field  stress  a  (t)  can  be  approximated: 


a(t) 


max 

”1 

a 


•  t 


0  £  t  s  t 

a 


a(t) 


t  <  t  S  t  , 
a  d 


The  derivative  differences  with  the  arguments  ^  are: 

t.  =0 

J 


£kf'  (0) 


max 


Af'  (t  )  =  a 

a  max 


^ax\t  t  I 
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(0  now  can 


t. 

J 


Af'(tci) 


<J 

max 


The  free-field  stress  as  a  function  of  the  circular  frequency 
be  expressed  by: 


a  (w) 


a  (to) 


Since  the  particle  velocity  is  bound  to  the  free-field  stress  by  the 
relationship: 


V 

where 

P 


a 


P 


C 

a 


=  density  of  the  free-field  medium 


C  =  longitudinal  seismic  velocity  of  the  free-field  medium 

or 

the  particle  velocity  as  a  function  of  the  circular  frequency  can  be  written; 


v(a)) 


max 


,  2 

2tr  o)  p 
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The  Fourier  transformation  delivers: 


v(t)  =  /  via,),  .““■‘-"‘'do, 
0 

or,  in  a  more  general  form: 


V  (t)  ~ 

0 

and 


v(t) 


iwt 

e  0(0 


DISPLACEMENT  POTENTIAL  AS  A  FUNCTION 
OF  THE  CIRCULAR  FREQUENCY 

The  transform  between  0(t)  and  ^  (co)  is: 


i  f  /  V  i(ax-(ot)  , 

0(t)  =  —  J(p(u))e  d(o 

0 

By  taking  the  partial  derivative  with  respect  to  x,  the  displacement 
equation  is: 

'i(t)  =  I  p(to)  ia  d<*) 

0 

By  taking  the  derivative  of  the  displacement  expression  with  respect  to 
time,  the  velocity  equation  is: 
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V(t)  = 

0 

and,  iwith  a  =  co/c  ,  the  form  i«; 

Qt 

/*»  1  7  /  »  2  i(ax-(ot)  j 

v(t)  =—  j^(a))<o  e'"  da) 

O  0 

By  comparing  this  equation  with  the  last  one  of  the  previous  section, 
the  frequency-dependent  potential  now  can  be  determined: 


0(60) 
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APPENDIX  F 


DERIVATION  OF  THE  STRESS  DISTRIBUTION  CAUSED  BY  A 
TRIANGULAR-SHAPED  PRESSURE  PULSE 


The  potentials  of  the  incident  triangular- shaped  pulse  that  propa¬ 
gates  in  the  x-  direction  are: 


(P 


00  i(ax  -  tot) 

(to)  e  dto 


=  0 


where 


p  (to )  =  Fourier  transform  of  0 
to  =  circular  frequency 


a 


(ti 


=  wave  number  of  the  longitudinal  wave. 


The  nonzero  longitudinal  potential  tp  can  be  expressed  by  a  Bessel 
function  series: 


w 

<p  =  j  <p  {(a) 


1  (ax- tot) 


00  00  ,  ^ 
dto  =  J  <p  (tolZ  jc^i’''^^(ar)cos(n0)e"^“jdto 


which  yields  stresses  as  follows: 


^  MM  2)  i",D^^^  +  A  +  B 

r(l)dy  1  J  “ J\  ®  “  1  nr  n  1  nr  n  1  nr  / 

®  n— UL' 

X  ^cos  n  0  e  dtojl 
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‘^0(l)dy  ' 
X 

^  r0(l)dy 
X 

®r(2)dy 

X 

®0(2)dy 

X 

^  re(2)dy 

X 


l.y^l  zf(0(.)c 

1  J  ^J\  o  n  1  nr  n  1  nr  n  1  nr  / 
CD  n=(\' 


(cos  n  0  e  do)) 


v-/ii 


0  (w)€  +A  +B 

’^o  n  1  nr  n  1  nr 


e(R)  j 

n  1  nr  / 


(sin  n  0 


-2 


W  UO 

/  2 

00  n=0 


M  +N 

n2  nr  n2 


+  R 

nr  n2  nr 


+  S 


5(R) 

n  2  nr  / 


\co8  n  0  e 


.2 


f;J 


M  +  N  +R  zE^^^+S 

n2  nr  n2  nr  n2  nr  n2  nr  j 


(  .  «  -iwt 

ysin  n  0  e 
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APPENDIX  G 


DERIVATION  OF  THE  ALLOWABLE  STRESSES  FOR  ROCK 


Figure  33.  Straight-line  envelope  of  Mohr's  circles. 


Notation; 

2 

K  =  cohesive  strength  in  kp/m 

0  =  angle  of  internal  friction 

or  j  =  major  principal  stress 

(j  ^  =  minor  principal  stress 

a  =  radial  stress 

r 

(7g  =  circumferential  stress 

=  shear  stress. 

ALLOWABLE  MAXIMUM  PRINCIPAL  STRESS 
IN  AN  ISOTROPIC  MEDIUM 

According  to  Mohr- Coulomb- Navier's  Theory,  the  envelope  of 
Mohr's  circles  is  the  straight  line; 


T 


=  tg  <0  •  a  +  K 


From  Figure  32,  it  can  be  read: 


T 

n 


cos  <p 


Q 

n 


a,  +02 


sin  tp 


sin  ifi  1^1  ^  ^2 
cos  ^  \  2 


n  (pj 


+  K 


Manipulations  lead  to  the  representation  of  the  envelope  related  to  the 
principal  stresses; 


a 


1 


2  K  cos  <p  ^  1  +  sin  p 
1  -  sin  (p  1  -  sin  (fi 


The  last  equation  means  that  for  a  certain  02  •  value  of  is  just 
the  allowable  limit.  Therefore,  the  equation  should  be  written: 


allowable 


2  K  cos  <p  ^  1  +  sin  <p 
1  -  sin  <p  1  -  sin  <p  ^2 


The  actual,  maximum,  principal  stress  has  to  yield  the  condition: 
®^1  ^allowable. 
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APPENDIX  H 


DERIVATION  OF  THE  SINGLE-SHOT  SURVIVABILITY 


Notation: 

Desired  ground  zero  is  the  origin  0  of  the  coordinate  system 

D  =  distance  between  desired  ground  zero  and  the  center  of  the 
structure  in  meters 

a  =  location  angle  for  the  center  of  the  structure  in  degrees 
L  =  length  of  the  tunnel  in  meters 

fi  =  angle  between  tunnel  direction  and  D  direction  in  degrees 
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R  =  failure  distance  of  the  tunnel  from  GZ 
a  =  standard  deviation  in  x  direction 

X 

a  =  standard  deviation  in  y  direction 

y 

P  =  probability  of  failure 

S  =  probability  of  survival  (survivability) 

The  idea  is  to  determine  the  target  area,  to  divide  the  target  area  in 
narrow  rectangles  of  width  Ax  and  of  length  y^  -  y^^  to  determine 
the  probability  of  failure  for  the  narrow  rectangles,  to  find  the  prob¬ 
ability  of  failure  over  the  entire  target  area  by  summation  of  the  prob¬ 
abilities  of  all  the  rectangles,  and  to  determine  the  survivability  out  of 
the  probability  of  failure. 

DETERMINATION  OF  THE  TARGET  AREA 


X  =  X  +  R  .  Bin  (a  -  /3) 
y^  =  -  R  .  cos  (a  -  /3) 
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Point  C^: 


V  -  +  R  •  ain  (a  -  B) 

Yq  =  -  R  •  cos  (a  -  P) 

Point  C^: 

X  =  X  -  R  •  sin  (a  -  p) 
^3  “2 

^2  ”2 

Point  C  . ! 

4 

-  R  .  sin  (a  -  P) 

4  1 

=  y„  +  R  •  cos  (a  -  p) 

C4  Hj 

BOUNDARIES 


Equation  of  circle  1 : 


Equation  of  circle  2: 


Equation  of  straight  line 


be  -Xc  j+Yc 


Equation  of  straight  line 
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PARTITION  OF  THE  TARGET  AREA  IN  NARROW  REaANGLES 


The  width  of  the  rectangles  is: 


The  divisor  A  determines  the  width  of  the  rectangles  in  relation  to  the 
weapon  accuracy  expressed  by  the  standard  deviation  in  the  x  -  direc¬ 
tion.  The  usual  problems  for  the  resultant  accuracy  (probability  result) 
was  found  to  be  sufficient  if  A  was  chosen: 

A  =  100 


gx 

100 


and  the  error  was  less  than  one  percent. 

To  get  the  length  of  the  rectangles  the  target-area  boundary  is 
divided  into  a  lower  boundary  BjC^C  and  an  upper  boundary 
B-C  C  B  (Figure  34).  The  length  of  rectangle  1  then  is: 


2 


) 


Because  of  the  different  functions,  the  y^  -  and  y^^  -  values  are  cal¬ 
culated  differently  over  the  x  -  range: 


(x„  -  R)  <  X  <  X. 


H 


y„=+[R^-(x.XH^)^ 


1/2 


s  X  i  X, 


y.,  = 


u  X  -  X 

^3  ^4 


(X  -  Xc  )  +  y^ 

4  4 
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The  cumulative  normal  distribution  function  is  used  to  represent 
the  single -shot  probability  in  military  targeting  (Figure  35).  Since  the 


Figure  35.  Normal  distribution  function. 
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integral  in  this  function  cannot  be  solved  by  the  direct  integration 
method,  numerical  values  were  obtained  by  series  approximations  and 
established  and  published  in  tables  in  References  11,  12,  and  13. 
Using  the  table  of  Reference  13,  four  series  functions  were  developed 
to  substitute  the  probability  function  in  the  range  from  z  =  .  00  to 
z  =  +00.  The  error  made  by  using  the  approximation  functions  is  less 
than  one  percent: 

1.95  s  z 


P(z)  =  1-0.  399 


0  <  z  <  1 . 95 


P(z)  =  0.841  +  0.242 
(-1.95)  <  z  s  0 

P(z)  =  0.  159  -  0.242 
z  s  (-1.95) 


P(z)  =  0.399  • 


The  independent  variable  z  is  measured  in  standard  deviations. 
In  targeting,  the  standard  deviation  is  expressed  in  a  length  unit.  To 
introduce  them,  the  length  values  x,  Ax,  y^,  and  y^  have  to  be  divided 
by  the  appropriate  standard  deviation: 


x 


X  +  Ax 

Z_  - - 

2  ox 
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^3  ay 


4  ay 

The  single -shot  probability  for  the  stripe  of  infinite  length  between  y 
and  y  is; 

Lj 


Py  =  P  ("3)  -  P  (%) 

Similarly,  the  single -shot  probability  for  the  stripe  of  infinite  length 
and  width  Ax  is: 

=  P  (z^)  -  P  (Zj) 

The  single -shot  probability  of  the  rectangle  x  •  (y  -  y  )  is: 

u  Xj 


p 

X 


Summation  of  the  whole  target  area  leads  to  the  single -shot  probability 
for  the  target: 


The  survivability  S  is  obtained  by: 


S  =  1  -  P 
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APPENDIX  I 


COST  AND  COST-EFFECTIVENESS  COMPUTATION 


access  tunnel  main  tunnel 


Figure  36.  Tunnel  section  for  example  problem. 


Notation; 

K 

C  , 
a2 

L 

Tu 


Tu 


Sh 


Tu 


Sh 


H 


=  cohesive  rock  strength,  in  kp/cm 

=  seismic  velocity  of  liner  material,  in  meters  per 
second 

=  length  of  main  tunnel,  in  meters 
=  length  of  access  tunnel,  in  meters 
=  inner  radius  of  access  tunnel,  in  meters 
=  inner  radius  of  tunnel,  in  meters 
=  inner  radius  of  access  shaft,  in  meters 
=  outer  radius  of  tunnel,  in  meters 
=  outer  radius  of  access  tunnel,  in  meters 
=  outer  radius  of  access  shaft,  in  moters 
=  depth  of  main  tunnel,  in  meters 
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] 


I 


>i 


w 


'Em 


Ex 


'Dr 


Co 


'St 


'StCy 


F 

C 


Tot 


'Tot  mach 


C 

C 

C 

C 


us 

Space 

Space  mach 

eff 

eff  mach 


percentage  of  longitudinal  reinforcement 

percentage  of  circumferential  reinforcement 

percentage  of  web  reinforcement 

cost  adjust  factor  for  labor  costs 

cost  adjust  factor  for  underground  equipment 

cost  adjust  factor  for  energy 

cost  adjust  factor  for  explosives 

cost  adjust  factor  for  drill  material 

cost  adjust  factor  for  concrete 

cost  adjust  factor  for  reinforcement  steel 

cost  adjust  factor  for  steel  liner 

cost  factor  for  shafts 

total  structure  costs  (dollars)  for  conventional  excava¬ 
tion  in  rock 

total  structure  costs  (dollars)  for  machine  excavation  in  rock 
usable  volume 

3 

structure  costs  per  usable  volume  in  rock  ($/m  ) 
for  conventional  excavation 

3 

structure  cost  per  usable  volume  in  rock  ($/m  ) 
for  machine  excavation 

cost  effectiveness  for  conventional  excavation 
cost  effectiveness  for  machine  excavation 
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The  cost  equations,  compiled  here,  are  based  upon  figures  in 
Reference  14. 

LABOR  COSTS 

The  labor  costs,  given  in  dollars  per  lineal  meter,  comprise 
costs  for  excavating  and  dumping  in  the  vicinity  of  the  tunnel  portal: 

Dry,  massive,  moderately  jointed  or  dry,  intact  rock: 

2 

K  i  200  kp/cm 

b  <  2. 30  m 

C.  =  E-  •  (88.  1  •  -  244.  1  •  b  +  450.  1) 

l-t  Lj 

b  2  2.30m 

C,  =  E,  (41.6  *  b^  -  91. 8  •  b  +  394.  1) 

Lu  Lj 

Dry,  stratified  or  schistose  rock: 

200  kp/cm^  >  K  i  100  kp/cm^ 
b  <  2. 30  m 

C-  =  E-  (55.4  •  b^  -  153.7  *  b  +  355.0) 

Lj  Lj 

b  2  2. 30  m 

C  =  E,  (43.3  •  b^  -  120.1  *  b  +  395.0) 

L  Xj 

Dry,  moderately  blocky  and  seamy  rock: 

2  2 
100  kp/cm  >  K  a  50  kp/cm 

b  <  2. 30  m 

C  =  E  •  (88.1  •  -  254.  4  •  b  +  438.0) 

Lj  Lj 

b  &  2.30m 

C  =  E-  •  (57.7-  b^  -  192.  8  *  b  +  504.0) 

Xj 
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Dry,  very  blocky  aiui  seamy  rock: 

2 

50  kp/cm  K 

b  <  2 . iO  m 

C,  =  E,  ■  (65.  5  •  b^  -  192.  7  •  b  t  416.0) 

L  ij 

h  >  2 .  3  0  ni 

=  E,  (44.  8  •  b'^  -  96.  3  •  b  4  364.  0)  . 

i-l 

UNDERGROUND  EQUIPMENT  COSTS 

Underground  equipment  includes: 

1.  Drill  jumbos,  drills,  mucking  machines ,  cars,  locomotives, 
compressors,  ventilator  fans,  smaH-miscellaneous  equip¬ 
ment,  and 

2.  Pipes  and  tracts. 

A  write-off  rate  of  15  percent  per  month  is  taken  into  account  for  the 
first  group.  This  rate  covers  depreciation  and  maintenance.  The 
monthly  costs  were  reduced  to  unit  costs  in  dollars  per  lineal  meter. 
The  second  group  of  equipment  items  could  be  related  directly  to  the 
tunnel  length. 

Dry,  massive,  moderately  jointed  or  dry  intact  rock: 

K  2  200  kp/cm^ 
b  <  1 . 8  m 

C_  =  E_  (13.6  •  b  +  144.0) 

£  £ 

1.8m<  b<  2.75m 

C_  =  E_  (32.2  ♦  b  +  195.  1) 

£  £ 

2.  75  m  s  b 

C_  =  E„  (64.4  •  b  +  160.8) 

£  £ 
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Dry,  stratified  or  schistose  rock; 

2  2 
200  kp/cm  >  K  >  100  kp/cm 

b  <  1 . 8  ni 

C  =E  •  154.1 

1.8m  ^  b<  2.75  m 

C_  =  E_  (32.2  •  b  +  165.  6) 

rL<  £ 

2.  75  m  ^  b 

C_  =  E_  (57.  3  •  b  +  141. 0) 

£  £ 

Dry,  moderately  blocky  and  seamy  rock; 

2  2 
100  kp/cm  >  K  s  50  kp/cm 

b  <  1.8m 

•  157.4 

E  E 

l.Sms  b<  2.75  m 

C_  =  E_  (32.2  •  b  +  175.5) 

£  £ 

2.  75  m  s  b 

C_  =  E_  (64.  5  •  b  +  134.  8) 

£  £ 

Dry,  very  blocky  and  seamy  rock; 

2 

50  kp/cm  >  K 

b  <  1.8m 

=  E^  •  164.0 
E  E 

1.8m^b<  2.75m 
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C  =  E  (48.  4  •  b  i  1 45. 0) 
E  E 


2.75  m  ^  b 


C  =  E^  (71.6  •  b  +  131.  1) 
E  E 


ENERGY  COSTS 

The  energy  coete  given  in  doliare  per  iineal  metera  are  baaed  on 
a  unit  price  of  1  cent  per  horsepower  hour; 

Dry,  massive,  moderately  jointed  or  dry,  intact  rock; 

.  2 

K  s;  200  kp/cm 

b  <  1.8m 

C  =  ■  14. 9 

En  En 

1.8ms  b  <  2.  75  m 

C  =  E^  (5.9  •  b  +  4.  3) 

^En  En 

2.  75  m  s  b 

C  =  E  (4.  6  ‘  b  +  10.  5) 

En  En 

Dry,  stratified  or  schistose  rock: 

2  2 
200  kp/cm  >  K  2  100  kp/cm 

b  <  1.8m 

=E^  •  13.1 
En  En 

1. 8  m  ^  b  <  2. 75  m 

=  E_  (3.  1  •  b  +  7.3) 

En  En 

2.75  m  ^  b 

C  =  E  •  (4. 1  •  b  +  6.  8) 

En  En  ^ 
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Dry,  moderately  blocky  and  seamy  rock: 

2  2 
100  kp/cm  >  k  s  50  kp/cm 

b  <  1.8m 


•  13.4 

En  En 


1.8m^  b<  2.75m 


=En=^En' 


2.  75  m  s  b 
'^En'^En- 

Dry,  very  blocky  and  seamy  rock: 
2 

50  kp/cm  >  K 
b  <  1.8  m 


=  E^  •  14. 1 
En  En 


1.8mS  b<  2.75m 


Sn  =  ^En  (^•8-b  +  5.9) 


2.  75  m  s  b 


COSTS  FOR  EXPLOSIVES 

The  costs  per  lineal  meter  are: 

Dry,  massive, moderately  jointed  or  dry  intact  rock 
K  i  200  kg/cm^ 
b  <  1. 0  m 
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=  0 
Ex 

l.Om^  b<  Z.Om 

•  (33.6  *6-21.  3) 

Ex  Ex 
2.0ms  b 

=  E^  (7.5  •  b^  -  3.4  •  b  +  22.2) 
Ex  Ex  ' 

Dry,  stratified  or  schistose  rock: 

2  2 
200  kp/cm  >  K  ^  100  kp/cm 

b  <  1.0m 

=  0 
Ex 

1.0msb<2.0m 

(25- ‘  -  14.  4) 

2.0  m  s  b 

•  (4.  1  *  b^  +  4.2  •  b  +  10.  5) 

Ex  Ex 

Dry,  moderately  blocky  and  seamy  rock: 

2  2 
100  kp/cm  >  Ks  50  kp/cm 

b  <  1.0  m 

=  0 
Ex 

1.0msb<2.0m 

<=Ex  =  =Ex'  P^-’-b->5.8) 


2.0  m  s  b 

'^Ex^^Ex'  <S-0-b'41.1  •  b4.5.1) 
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Dry,  very  blocky  and  seamy  rock: 

2 

50  kp/cm  >  K 

b  <  1.0m 

=  0 
Ex 

1.0m^b<2.0m 

=  (2‘-5-b-U.4) 

2.0  m  <  b 

=  (3.0-  bS  5.9  -5  +  7.6)  . 

COST  FOR  DRILL  BITS  AND  RODS 

The  costs  for  drill  bits  and  rods  are  as  follows: 
carbide  inset  bits— $14. 20  each 
drill  rods— $15.00  each 

The  unit  costs  are  dimensioned  in  dollars  per  lineal  meters.  They  are 
for  all  kinds  of  rocks: 

*  (3.  3  •  b  -  1.2)  . 

Dr  Dr  ' 


EXCAVATION  COSTS 

The  sum  of  the  labor,  equipment,  energy,  explosive,  and  drill 
material  costs  establishes  a  subtotal  of  the  excavation  costs: 


Sub 


=  C,  +  +  C_  +  C_  +  C 


En 


Ex 


Dr 


For  engineering,  overhead,  and  contingencies,  25  percent  of  the  sub¬ 
total  is  added.  For  profits,  15  percent  of  the  subtotal  is  added.  In 
addition  to  these  costs,  expenses  for  miscellaneous  items  are  to  be 
taken  into  account.  So,  the  excavation  costs  in  dollars  per  lineal  meter 
are: 
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I 


Li 


i 


EC  Sub  M 

where  C. .  is  the  miscellaneous  cost; 
M 


LIMfR  COSTS 

The  tunnels  might  be  lined  in  plain  concrete,  reinforced  concrete, 
or  steel.  The  seismic  velocity  of  the  liner  is  distinguishing  in  the  pro¬ 
gram  when  concrete  or  steel  is  to  be  considered  in  the  cost- computation 
subroutine.  If  the  seismic  velocity  liner  equals  or  is  less 

than  4,  500  meters  per  second,  control  is  given  to  the  calculation  of  the 
costs  of  reinforced  concrete  liner.  If  the  seismic  velocity  exceeds 
4,  500  meters  per  second,  the  costs  for  a  steel  liner  are  computed. 


REINFORCED  CONCRETE  LINER 

The  amount  of  reinforcement  steel  is; 


St 


=  W  [("<  ^  Pe)  +  2  •  “  •  p*  •  (» -  »)] 


The  amount  of  concrete  is; 
A 


Co 


=  3.  14  [(o.  15  +  b)^  -  a^j  -  Ag^ 


The  costs  for  the  concrete  are  based  on  a  unit  price  of  $40.  00  per 
cubic  yard  in  place.  The  costs  for  the  reinforcement  steel  are  based 
on: 


$0.  30  per  pound 


1 

'Co 

=  • 

Co 

A  • 

Co  0. 76 

St 

=  E  • 
St 

A_.  •  7.  85  •  2204.  6  •  0.  30 
ot 

99 


The  costs  per  lineal  meter  for  a  reinforced  concrete  liner  are: 


Li 


=  C  +  C 
Co  St 


If  plain  concrete  is  used,  the  costs  are  obtained  by  setting: 

=  Pe  =  Pw  =  “ 


The  costs  for  steel  liner  consist  of  the  costs  for  steel  and  Gunite 
between  rock  and  steel  liner.  The  basic  unit  prices  are: 

$0.44  per  pound  for  steel  in  place 
$3.  95  per  cubic  foot  for  Gunite  in  place 


Li 


=  E__  (b^  -  a^)  •  3.  14  .  2204.  6  •  0.  44 


+  •  6.  28  •  b  •  0.05  •  35  •  3.  95  . 

Co 


TOTAL  STRUCTURE  COSTS 


The  sum  of  the  excavation  costs  and  the  liner  costs  furnishes  the 
total  structure  costs  per  lineal  meter: 


'Str 


Li 


The  protective  underground  structure  consists  usually  of  the  main 
tunnel,  access  tunnels,  and  access  shafts.  The  costs  of  the  access 
structures  can  be  calculated  with  the  previously  given  cost  equations, 
as  well  as  the  costs  of  the  main  structure.  For  the  access  structures, 
the  pertaining  radii  are  to  be  introduced: 


Access  tunnel: 


'n,  =  C-.  (a  =  a  ,  b  =  b  ) 

Tu  Str  Tu  Tu 


Access  shaft: 


C  =C  (a  =  a  ,b=b  )*F 

^Sh  ^Str  '  Sh  ’  Sh' 
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The  factor  F  in  the  cost  equation  for  the  shaft  brings  into  account 
that  a  shaft  is  more  expensive  than  a  tunnel  of  the  same  diameter.  In 
general,  F  =  1.  5  .  The  total  costs  in  dollars  for  the  access  structures 
are: 


AcTu 


^AcSh 


=  (H  +  b  +  3)  • 


The  access  costs  are; 


Ac 


^AcTu  ^  ^AcSh 


The  total  costs  for  the  entire  protective  structure  are: 
C. 


Tot 


=  C„  •  L  +  C  .  . 

Str  Ac 


SPECIFIC  STRUCTURE  COSTS 

In  many  application  cases,  an  estimation  of  the  total  structure 
cost  is  insufficient.  Very  often  specific  costs  must  be  made  available 
for  economic  comparisons— cost  related  to  usable  volume; 

V  =  L  •  3.  14  •  a^ 

us 

is; 

Cc 

Space 

The  effectiveness  and  the  cost  of  a  structure  are  measured  in  the  ex¬ 
pression  of  the  cost  effectiveness: 


Tot 


us 
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APPENDIX  J 


DESCRIPTION,  FLOW  CHART,  AND 
LISTING  OF  THE  MAIN  PROGRAM 


A  flow  chart  for  the  computer  program  is  given  in  Figure  37. 
The  arrangement  of  the  program  can  be  understood  and  followed  by 
using  the  flag  "NPAS"  as  a  guide.  Prior  to  setting  "NPAS"  to  one, 
the  Input  is  read  from  five  cards; 


Card 

1:  Format:  2F10.  0,  2F5.  0,  SE  10.6 

Cols 

1-10 

Specif,  gravity  of  media 

(no  units) 

Cols 

11-20 

Specif,  gravity  of  liner 

(no  units) 

Cols 

21-25 

Poisson's  ratio  of  media 

(no  units) 

Cols 

26-30 

Poisson's  ratio  of  liner 

(no  units) 

Cols 

31-40 

Elastic  modulus  of  media 

(kp/m^) 

Cols 

41-50 

Elastic  modulus  of  liner 

(kp/m^) 

Cols 

51-60 

Radius  of  inside  liner 

(meters) 

Cols 

61-70 

Radius  of  outside  liner 

(meters) 

Cols 

71-80 

Depth  of  cylinder  center 

(meters) 

Card  2:  Format:  7E10.  6,  F5.  0,  15 

Cols 

1-10 

Tangent  of  phi  of  media 

(no  units) 

Cols 

11-20 

Cohesive  strength  of  media 

(kp/m^) 

Cols 

21-30 

Tensile  strength  of  media 

(kp/m^) 

Cols 

31-40 

Constant  "L"  of  liner 

(kp/m^) 

1  2  1 

Cols 

41-50 

Constant  "M"  of  liner 

(kE)  .  /IE_) 

'2'  '  kp  ^ 

m  ^ 

Cols 

51-60 

Constant  "N"  of  liner 

(no  units) 

Cols 

61-70 

Tensile  strength  of  liner 

(kg/m^) 
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c 

r 

' 

Cols 

71-75 

Nominal  Weapon  Yield 

(KT) 

Cols 

76-80 

Number  of  weapon  type 

(no  units)  i 

V 

Card 

3:  Format:  8F10.  0 

! 

Cols 

1-10 

Distance  from  cylinder  to  target 

(meters) 

Cols 

11-20 

Length  of  cylinder 

(meters) 

Cols 

21-30 

Angle  alpha 

(degree) 

Cols 

31-40 

Angle  beta 

(degree) 

1 

Cols 

41-50 

Standard  deviation  in  "X" 

(meters) 

Cols 

51-60 

Standard  deviation  in  "Y" 

(meters) 

Cols 

61-70 

Effect  longitudinal  seismic  velocity 
between  GZ  and  structure 

(meters 
per  second) 

Cols 

71-80 

Adjustment  factor  for  labor  cost 

(no  units) 

Alpha: 

Location  of  cylinder  from  target 

Beta: 

Orientation  of  cylinder.  . .  direction 

Standard  Deviation;  Range  at  which  40%  of  shots  fall  inside. 

Card 

4:  Format:  8F10. 0 

Cols 

1-10 

Radius  inside  liner  of  access  tunnel 

(meters) 

Cols 

11-20 

Radius  outside  liner  of  access 
tunnel 

(meters) 

Cols 

21-30 

Access  tunnel  length 

(meters) 

Cols 

31-40 

Radius  inside  liner  of  access  shaft 

(meters) 

Cols 

41-50 

Outside  liner  of  access  shaft 

(meters) 

Cols 

51-60 

Ratio  of  longitudinal  reinforcement 

(percent) 

1 

Cols 

61-70 

Ratio  of  circumferential  reinforce¬ 
ment 

(percent) 

1 

? 

!►  1 
i 

Cols 

71-80 

Ratio  of  web  reinforcement 

(percent) 

» 

i.  J 

*  Number  of  weapon  type;  "1"  Contact  surface  burst,  "2"  true  surface 

burst,  "3"  fully  buried  detonation  (nominal  weapons  yield  =  effective 

weapons  yield).  ! 

1 

i 
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Card 

5:  Format 

8F10.  0 

Cols 

1-  10 

Adjustment  factor  for  equipment 
cost 

(no  units) 

Cols 

11-20 

Adjustment  factor  for  power 

(no  units) 

Cols 

21-30 

Adjustment  factor  for  explosives 

(no  units) 

Cols 

31-40 

Adjustment  factor  for  drill  bits 
and  rods 

(no  units) 

Cols 

41-50 

Adjustment  factor  for  concrete 

(no  units) 

Cols 

51-60 

Adjustment  factor  for  reinforce¬ 
ment  steel 

(no  units) 

Cols 

61-70 

Miscellaneous 

(no  units) 

Cols 

71-80 

Adjustment  factor  for  steel  liner 

(no  units) 

.NPAS  =  1. 

The  static  stresses,  the  principal  stresses,  and  the  yield  stresses 

are  computed  for  inside  the  liner.  A  yield  check  is  made.  If  yielding 
does  not  occur,  then  it  is: 

NPAS  =  2. 


The  procedure  under  NPAS  =  1  will  be  repeated  for  the  interface 
for  both  liner  and  rock.  If  no  yielding  exists,  the  weapon  values  will  be 
computed.  The  horizontal  distance,  incremented  by  one  diameter,  is 
set  to  zero.  The  peak  free-field  dynamic  stress  is  determined  and 
checked  against  a  preliminary  yield: 


preliminary  = 


2K  •  cos 
1  -  sin  (p 


where,  in  accordance  with  the  notation  of  Appendix  G, 

2 

K  =  cohesive  strength  of  rock  in  kp/m 

<P  =  angle  of  internal  friction  in  degrees 

If  yielding  occurs,  the  horizontal  distance  will  be  increased  as  long  as 
the  check  against  the  preliminary  yield  is  negative.  The  application  of 
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the  preliminary  check  shortens  considerably  the  run  time  since  the 
dynamic  stresses  and  the  appropriate  strength  values  need  not  be 
computed  for  each  new  horizontal  distance. 

If  the  preliminary  yield  check  is  negative,  then  it  will  be: 

NPAS  =  3. 


The  pressure  pulse  is  computed.  The  frequency  values  for  the 
integration  are  established: 

NPAS  =  4, 

Static  and  dynamic  stresses  and  the  pertaining  yield  values  are 
computed  for  rock  and  liner  at  the  interface.  If  either  rock  or  liner 
yields,  the  horizontal  distance  is  increased  and  the  dynamic  part  is 
resolved: 

NPAS  =  5. 

The  procedure  of  NPAS  =  4  is  repeated  for  the  inside  of  the  liner. 
.If  it  yields  no  longer,  then  it  will  be: 

NPAS  =  6. 

Where  the  failure  occurred  first,  in  the  interface  or  at  the  inside 
of  the  liner,  is  checked.  The  horizontal  distance  is  moved  into  the  safe 
area  by  increasing  it  by  one  diameter.  The  survivability  and  the  eco¬ 
nomic  features  are  computed.  Then  a  sixth  card,  containing  information 
about  the  type  of  the  next  problem,  is  read  in: 

Card  6:  Format:  II. 

Only  one  value  occurs  in  Col.  1; 

"0"  All  data  have  been  run  and  the  program  calls  an  exit 
and  prints  "Job  Completed.  " 

"1"  A  new  problem  is  to  be  run  in  which  new  data  cards 
(1-6)  will  be  read. 

"2"  Only  new  data  cards  (3-6)  will  be  read.  This  allows  a 
change  in  the  cost  of  probability  programs:  one  change 
in  the  length  of  the  cylinder  or  distance  to  the  target. 

The  output  data  are  the  static  stresses  for  the: 

1.  Inside  liner. 


105 


2.  Interface  liner,  and 

3.  Rock  at  the  interface  for  13  angles  from  0  to  180  degrees. 
Other  output  data  include: 

1.  Total  stresses  for  13  angles  and  10  time  steps, 

2.  Minimum  survival  distance, 

3.  Probability  of  failure, 

4.  Survivability, 

5.  Total  cost  for  conventional  tunneling, 

6.  Total  cost  for  machine  tunneling,  and 

7.  Cost  per  usable  volume  per  unit  increment  of  survivability. 
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Figure  37— continued. 
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Figure  37— continued. 
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Figure  37— continued. 
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(U)  The  objective  of  this  study  wis  to  define  the  limit  survival  distances  from  ground 
aero  for  deep  underground  protective  structures  in  rock  and  to  formalize  a  methodology  for 
defining  the  cost  of  such  systems.  To  achieve  this  objective,  a  computer  program  was  developed 
and  subsequently  was  exercised  for  representative  tunnels  with  concrete  and  steel  liners  in 
sandstone  or  (panite. 

Limit  distance  contours  were  obtained  for  a  range  of  effective  seismic  velocities.  They 
show  that  there  are  optimum  depths  which,  if  exceeded,  will  increase  the  survival  distance  and 
dacraasa  the  survivability.  The  optimum  site  profile  consists  of  a  thick  soil  layer  over  a  strong 
basement  rock. 

Steel  liners  provide  greater  survivability  and  lower  cost  than  concrete  liners. 
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yMd  weapons  than  for  larger  ones.  It  is  apparent  from  the  resuitt  that  deep  underground 
protective  structures  can  be  designed  for  high  degrees  of  survivability. 
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